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Abstract 

We  examine  two  problems  concerning  semilinear  elliptic  equations.  We  consider 
single  equations  of  the  form  An  =  p(x)ua  +q(x)ufi  for  0  <«</?<  1  and  systems 
A  u  =  p(\x\)f  (v),  Av  =  q(\  x  \  )g(u)  ,  both  in  Euclidean  n  -space,  n  >  3  .  These  types  of 
problems  arise  in  steady  state  diffusion,  the  electric  potential  of  some  bodies,  subsonic 
motion  of  gases,  and  control  theory.  For  the  single  equation  case,  we  present  sufficient 
conditions  on  p  and  q  to  guarantee  existence  of  nonnegative  bounded  solutions  on  the 
entire  space.  We  also  give  alternative  conditions  that  ensure  existence  of  nonnegative 
radial  solutions  blowing  up  at  infinity.  Similarly,  for  systems,  we  provide  conditions  on 
p,q,f,  and  g  that  guarantee  existence  of  nonnegative  solutions  on  the  entire  space.  The 
main  requirement  for  /  and  g  will  be  closely  related  to  a  growth  requirement  known  as 
the  Keller-Osserman  condition.  Further,  we  demonstrate  the  existence  of  solutions 
blowing  up  at  infinity  and  describe  a  set  of  initial  conditions  that  would  generate  such 
solutions.  Lastly,  we  examine  several  specific  examples  numerically  to  graphically 
demonstrate  the  results  of  our  analysis. 
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ENTIRE  BLOW-UP  SOLUTIONS  OF  SEMILINE AR  ELLIPTIC  EQUATIONS  AND 


SYSTEMS 


I.  Introduction 

Our  work  is  divided  into  two  parts.  First,  we  examine  the  semilinear  elliptic  equation 

A u  =  p(x)ua  +  q(x)up ,  x  e  R"  (1-1) 

for  n  >  3  .  This  problem  was  the  focus  of  the  thesis  by  Smith  [26]  in  which  multiple  existence 
and  nonexistence  results  were  obtained.  We  examine  this  problem  only  briefly  in  an  effort  to 
close  a  few  remaining  gaps.  The  arguments  we  present  are  similar  to  those  used  in  [26]  and  the 
works  referenced  therein. 

The  second  part  of  our  study  is  dedicated  to  semilinear  elliptic  systems.  Namely,  we 
examine 


| Am  =  p(\  x  |)/(v), 
1  Av  =  q{ |  x  |)g(u), 


where  again  n>  3  .  Here  we  are  restricting  our  problem  to  the  radial  case.  Several  authors  have 
studied  this  system  for  mono  tonic  /  and  g  .  We  will  consider  a  more  general  system  where  / 
and  g  may  be  non-mono  tonic.  To  our  knowledge,  no  results  exist  for  system  (1.2)  when  /  and 
g  are  not  monotone.  This  generalization  presents  several  unique  challenges  that  we  must 
address. 
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For  both  single  equations  and  systems,  we  are  primarily  concerned  with  proving  the 
existence  of  large  solutions.  A  solution  to  ( 1 . 1 )  is  large  if  u(x)  — »  oo  as  x  — >■  dfl .  If  Q  =  R" ,  we 
require  u(x)  — >  oo  as  |  x  |— >  oo  .  The  latter  case  is  called  an  entire  large  solution.  Similarly,  a 
large  solution  to  system  (1.2)  is  a  solution  (w,v)  such  that  u(  \  x  |)  — >  oo  and  v(|  x  |)  — »  oo  as 
|  x  |— >  dQ. ,  and  an  entire  large  solution  is  one  such  that  u(\  x  |)  — »  oo  and  v(|  x  |)  — >  oo  as  |  x  |— >  oo  . 

The  arguments  developed  in  our  analysis  and  the  mathematical  foundations  of  elliptic 
theory  may  be  applied  to  multiple  problems  in  a  wide  variety  of  technical  fields.  Elliptic 
equations  similar  to  those  we  analyze  here  are  related  to  steady-state  reaction-diffusion,  subsonic 
fluid  flows,  electric  potentials  of  some  bodies,  and  control  theory. 

For  example,  [17]  describes  a  general  stochastic  diffusion  process  with  feedback  controls. 
The  controls  are  to  be  designed  so  that  the  state  of  the  system  is  constrained  to  some  region. 
Finding  optimal  controls  is  then  shown  to  be  equivalent  to  finding  large  solutions  for  a  second 
order  semilinear  elliptic  equation.  As  another  example,  [8]  models  the  steady  state  of  non-linear 
heat  conduction  through  a  2-component  mixture  with  a  system  of  semilinear  elliptic  equations 
similar  to  the  systems  we  study  here 

These  are  just  two  examples  of  application.  The  mathematical  methods  and  ideas  we  use 
and  develop  here  are  applicable  in  many  scientific  and  engineering  disciplines.  Before  we 
discuss  preliminaries  for  our  work,  let  us  first  examine  results  obtained  by  others  in  this  field. 
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II.  Background 


Bieberbach  [3],  in  1916,  was  the  first  to  study  large  solutions  to  the  semilinear  elliptic 
problem 


A u  =  f(u),  reQ  (2.1) 

where  / (u)  =  e“ .  Since  that  time,  many  authors  have  studied  related  problems  for  single 
equations  and  systems.  In  this  section  we  will  present  only  the  most  recent  and/or  most  relevant 
accomplishments  which  have  led  to  our  study. 

2.1  Single  Equations 

In  1957,  Keller  [11]  and  Osserman  [21]  established  necessary  and  sufficient  conditions 
for  the  existence  of  solutions  to  (2.1)  on  bounded  domains  in  M" .  They  showed  that  large 
solutions  exist  on  Q  if  and  only  if  /  satisfies 


dt  <  oo. 


(2.2) 


This  requirement  is  sometimes  referred  to  as  the  Keller-Osserman  condition  and  continues  to  be 
significant  in  current  studies.  Indeed,  we  will  require  several  of  our  functions  to  satisfy  this 
Keller-Osserman  condition. 

Bandle  and  Marcus  [2]  later  examined  the  equation 

A  u  =  p(x)f(u)  (2.3) 
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for  /  non-decreasing  on  [0,  oo)  and  proved  the  existence  of  positive  large  solutions  provided  the 


function  /  satisfies  (2.2)  and  p  is  continuous  and  strictly  positive  on  Q  .  Lair  [12]  showed  the 
same  results  hold  for  (2.3)  when  p  is  allowed  to  vanish  on  large  parts  of  Q ,  including  its 
boundary.  Lair,  Proano,  and  Wood  [16]  relaxed  the  monotonicity  condition  on  /  by  requiring 
that  there  exist  some  nonnegative,  nondecreasing,  Holder  continuous  function  g  and  positive 
constants  M  and  s0  such  that 

£•(£)<  f(s)<Mg(s )  for  all  s>s0.  (2.4) 

We  note  that  our  work  on  systems  is  a  similar  generalization;  however,  the  arguments  required 
for  our  analysis  are  very  different  from  those  used  in  the  single  equation  case. 

Many  authors  have  examined  more  specific  forms  of  (2.3).  The  equation 

Au  -  p(x)ur  (2.5) 

has  been  of  particular  interest.  Cheng  and  Ni  [4]  considered  the  superlinear  case  (y  >  1)  and 
proved  that  (2.5)  has  large  solutions  on  bounded  domains  provided  p  is  strictly  positive  on  <9Q  . 
Lair  and  Wood  [14]  generalized  this  to  allow  p  to  vanish  on  large  portions  of  Q  including  its 
boundary.  They  also  showed  the  existence  of  an  entire  large  solution  to  (2.5)  provided  that 

oo 

[  r  max  p(x)dr  <  oo  .  (2.6) 

J  |x|=r  ' 

0  1  1 

This  is  a  weaker  condition  compared  to  the  requirements  in  [4], 
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Fewer  results  are  known  for  the  sublinear  case  ( 0  <  y  <  1 )  of  (2.5).  In  [15],  Lair  and 
Wood  proved  that  entire  large  radial  solutions  of  (2.5)  exist  if  and  only  if 

oo 

J  rp(r)dr  =  oo  .  (2.7) 

0 

They  also  demonstrated  that  for  a  bounded  domain  Q ,  (2.5)  has  no  positive  large  solution  when 
p  is  continuous  in  Q  .  In  addition  to  large  solutions,  Lair  and  Wood  considered  entire  bounded 
solutions.  They  proved  entire  bounded  solutions  of  (2.5)  exist  when  (2.7)  holds  and  p  is  locally 
Holder  continuous.  Also,  nonnegative,  entire  bounded  solutions  do  not  exist  for  (2.5)  when 

oo 

[  r  min  p(x)dr  =  oo  .  (2.8) 

J  \x\=r 

0  1  1 

Smith  [26]  then  considered  (1.1),  which  is  a  multi-term  adaptation  of  the  single  term 
equation  (2.5).  For  the  superlinear  (1  <  a  <  J3)  and  mixed  (0  <  a  <  1  <  J3)  cases,  results 
comparable  to  those  for  single  equations  in  [14]  were  obtained.  In  the  sublinear  case,  Smith  [26] 
proved  the  existence  of  entire  large  radial  solutions  and  showed  the  existence  of  entire  bounded 
solutions  for  the  nonradial  problem.  However,  both  of  these  proofs  require  0  <  a  <  ft  <  1 .  The 
proofs  do  not  hold  when  fd  - 1 .  Table  2-1  is  a  summary  of  the  results  from  [26].  Our  work  on 
single  equations  removes  condition  (f)  from  the  table,  thereby  closing  the  gap  for  [5  —  1 . 
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Table  2-1.  Existence  of  Solutions  for  Multiple  Term  Single  Equation  (1.1) 


Superlinear/Mixed 

0  <a<  P,  P>  1 

Sublinear 

0  <  a  <  (3  <  1 

Sublinear 

0  <  a  <  P  <  1 

Bounded 

Domain 

Qc=M” 

Entire 

Domain 

Q  =  M" 

Bounded 

Domain 

Qcl" 

Entire 

Domain 

Q  =  M" 

Entire 

Domain 

Q  =  M'! 

Large 

Solutions/ 

Requirements 

Yes 

a 

Yes 

a,b 

No 

Yes 

c,d,f 

Bounded 

Solutions 

Yes 

b,f 

No 

e 

Au  =  p(x)ua  +  q(x)ut 


a)  p,  q  are  c-positive  (see  Definition  4-2) 

oo  oo 

b)  |  r  max  p(x)dr  <  oo  and  [  r  max  q(x)dr  <  co 

J  bc=r  *  x=r 

0  1  1  0  1  1 

c)  radial  case  only:  p{x)  =  /?(|x|)  and  q{x)  =  ^(|x|) 

oo  oo 

d)  [  r  max  p(x)dr  =  co  or  \r  max  q(x)dr 

J  |jcl=r  J  Ul=r 


e)  [  r  min  p(x)dr  =  co  or  [  r  min  q(x)dr 

J  |x|=r  J  |x|=r 


=  00 


=  00 


f) 


2.2  Systems 

Next  we  shall  discuss  the  background  of  semilinear  elliptic  systems  of  the  type  given  in 
(1.2).  While  systems  are  a  natural  extension  of  single  equations  and  occur  in  many  of  the  same 
areas  of  application,  not  all  the  methods  employed  to  study  them  carry  over.  For  instance,  we  do 
not  have  a  meaningful  maximum  principle  for  systems.  Additionally,  the  barrier  method,  which 
we  present  later,  is  a  common  method  for  proving  the  existence  of  solutions  for  single  equations. 
However,  a  similar  method  for  systems  exists  only  when  rather  restrictive  conditions  are  placed 
on  the  growths  of  the  given  functions.  For  example,  if  we  consider  systems  built  from  equations 
comparable  to  (2.1),  we  have 


|A  u  =  f(u,v), 
[Av  =  g(u,v). 
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A  barrier  method  similar  to  that  for  single  equations  exists  only  when  /  is  monotonic  in  v ,  and 
g  is  mono  tonic  in  u  .  Still,  this  area  has  been  well  studied  by  several  authors.  Most  related  to 
our  study,  are  the  works  of  Lair  and  Wood  [13],  Cirstea  and  Radelescu  [6],  and  Peng  and  Song 
[22]. 

Lair  and  Wood  [13]  examined  the  system 


| Am  =  p(\x\)va, 
[Av  =  q(\x\)u/1, 


xeQci?" 


(2.9) 


where  n>  3  .  They  proved  the  existence  of  entire  nonnegative  solutions  and  characterized  the  set 
of  central  values,  ScR+xl+,  for  these  solutions.  By  central  values,  we  mean  (m(0),  v(0)) 
where  u,v  are  solutions  of  (2.9).  Lair  and  Wood  proved  S  is  closed,  bounded,  and  convex. 

They  also  further  geometrically  characterized  this  set  by  describing  bounds  for  S  .  Finally,  they 
proved  that  large  solutions  exist  for  central  values  that  lie  in  the  closure  of  {( a,b )  f|  dS :  a,b  ^  0} . 
We  note  that  for  a,  J3  >  1 ,  they  required 


jtp(t)<  co,  j"  tq(t)  <  oo  .  (2.10) 

0  0 

Both  Peng  and  Song  [22]  and  Cirstea  and  Radulescu  [6]  generalized  their  work  to  show 
comparable  results  for  system  (1.2)  where  /,g  e  C[0,oo)  are  nonnegative  monotonic  functions. 
The  approaches  taken  in  each  paper  were  slightly  different,  and  each  pair  of  authors  required 
different  additional  constraints  on  /  and  g  .  In  [6]  several  cases  were  considered.  When 
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limlWM  =  0forallc>0,  <2U) 

f->a>  f 

system  (1.2)  has  entire  solutions,  all  of  which  are  bounded  when  (2.10)  holds.  Further,  entire 
solutions  exist  and  are  large  when  neither  inequality  in  (2.10)  hold.  Both  papers  then  considered 
(1.2)  where  the  nonnegative  functions  p,qeC[  0,oo)  satisfy  (2.10),  and  functions  f,geC[  0,oo) 
are  nondecreasing,  satisfy  the  Keller-0 ssennan  condition  (2.2),  and 

/( 0)  =  g( 0)  =  0 ,  f(s )  >  0,  gfy)  >  0  for  5  >  0 .  (2. 12) 

In  [22],  these  functions  were  further  required  to  be  convex,  while  in  [6],  there  were  additional 
conditions  /,geC'[ 0,oo)  and 


liminf  =  a  >  0  .  (2  131 

—  g(s)  1  ' 

In  these  cases,  both  pairs  of  authors  showed  existence  of  entire  large  solutions  and  characterized 
the  set  of  central  values  for  which  such  solutions  exist. 

Our  results  on  systems  further  generalize  this  problem  by  allowing  /  and  g  to  be  non¬ 
monotonic.  Instead,  we  will  require  the  function  G  ,  given  by 

G(s)  =  mini  min  f(t),  min  g(t)i,  0<s,  (2.14) 

l.  s<t  s<t  ) 

to  satisfy  the  Keller-Osserman  condition  in  (2.2).  Note  this  automatically  implies  /  and  g  must 
satisfy  this  condition  as  well. 
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Now  we  discuss  preliminary  theory  that  will  be  necessary  for  both  single  equations  and 
systems. 

2.3  Preliminaries 

For  readability,  we  shall  reserve  the  terms  lemma,  theorem,  and  corollary  for  results 
which  we  justify  in  our  work.  Any  result  taken  directly  from  another  reference  shall  be  labeled 
as  Preliminary  X-X. 

The  Arzela-Ascoli  Theorem  is  our  main  tool  for  examining  systems;  however,  it  is  also 
needed  for  a  single  equation  result.  Here  we  give  the  Arzela-Ascoli  Theorem  and  several 
necessary  definitions. 

Definition  2-1  A  subset  K  of  a  nonned  space  X  is  called  precompact  (sometimes  called 
relatively  compact  or  conditionally  compact)  if  its  closure  K  is  compact  in  the  norm  topology  of 
X. 


Preliminary  2-1  (Theorem  1.34  of  [1])  (Arzela-Ascoli  Theorem)  Let  Q  be  a  bounded  domain  in 
M" .  A  subset  K  of  C(Q)  is  precompact  in  C(Q)  if 

(i)  There  exists  M  >  0  such  that  \<fi(x)  \  <  M  for  every  (/)  e  K  and  x  e  Q  .  That  is,  K 
is  uniformly  bounded. 

(ii)  For  every  £  >  0 ,  there  exists  8  >  0  such  that  \<f>(x)  -  <f>(y)\  <£  for  all  (j)  e  K , 
rjefi,  and  \x-  y\<  8 .  That  is,  K  is  equicontinuous. 
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To  apply  the  Arzela-Ascoli  Theorem,  we  will  often  show  that  a  sequence  of  functions  is 
both  unifonnly  bounded  and  equicontinuous.  We  will  need  to  accomplish  this  for  several 
different  sequences  throughout  our  arguments.  However,  these  sequences  will  have  similarities, 
and  creating  the  following  lemma  will  greatly  streamline  our  work  in  later  sections.  This  lemma 
will  allow  us  to  conclude  a  particular  sequence  is  both  uniformly  bounded  and  equicontinuous 
only  by  showing  unifonn  boundedness. 

Lemma  2-1  Let  {uk }  be  a  sequence  of  functions  of  the  form 

r  t 

uk(r)  =  ak  + 1  tl~n^sn~lp{s)f  (vk(s))dsdt,  re[0,l?].  (2.15) 

0  0 

where  ak  e  1R+  ,  p,f  e  C[0,oo)  are  nonnegative,  and  {v’J  is  an  arbitrary  sequence  of 
nonnegative  continuous  functions  on  [0,f?] .  If  the  sequence  { vy. }  is  uniformly  bounded  on 
[0 ,R] ,  then  \uk }  is  equicontinuous  on  [0,  R\ . 

Proof.  Since  { vk  }  is  uniformly  bounded,  there  exists  M  such  that  vk  (r)  <M  for  all  k  and  all 
r  e  [0,  f?] .  Then  we  have 
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u'k{r)  =  r |  sn~1p(s  )f(vk  (s  ))ds 
0 

<r~nrn~l^p{s)f{vk{s))ds 

0 

r 

<  f  p(s)  max( f  (t))ds 

J  0  <t<M 

0 

R 

<  max(/(0)J  p(s)ds 

o 

S  C  <  00. 


Also,  we  clearly  have  0  <  u[{r )  .  Thus,  {uk}  is  equicontinuous  on  [0  ,R] ,  and  our  proof  is 
complete.  □ 


These  are  the  basic  ideas  needed  for  both  the  single  equation  and  system  cases.  Now  we 
turn  our  attention  to  preliminaries  required  exclusively  for  our  single  equation  arguments. 
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III.  Single  Equations 


3.1  Preliminaries 

We  begin  by  presenting  a  result  on  the  barrier  method  or  upper- lower  solution  approach. 
This  method  is  well-known  for  equations  on  bounded  regions  (see  Theorem  2.3.1  of  Sattinger 
[24]).  However,  we  wish  to  apply  the  barrier  method  on  the  unbounded  domain  M" .  We 
therefore  use  the  following  result  (Preliminary  3-1)  from  Shaker  [25],  We  also  provide 
definitions  related  to  this  result. 

Definition  3-1  Let  /  be  a  function  defined  on  an  open  set  QcS".  For  0 <  a  <  1 ,  we  say  that 
/  is  Holder  continuous  with  exponent  a  ,  written  /  e  Ca  on  Q  ,  if  for  every  x,y  e  Cl,  there 
exists  a  nonnegative  constant  C  such  that 

\f(x)-f(y)\<C\x-y\a. 

When  a  =  1 ,  we  say  /  is  Lipschitz  continuous.  Further,  we  say  that  /  is  locally  Holder 
(Lipschitz)  continuous  on  QcS"  when  for  every  jeQ,  there  exists  an  open  ball  B(x,  r)  such 
that  /  restricted  to  B{x,  r)  is  Holder  (Lipschitz)  continuous. 

Preliminary  3-1  (Lemma  3  of  [25])  (Lemma  on  Barrier  Method)  Let  ux,u2 :  M”  — >  R , 
ux(x)  >  u2{x)  for  all  x  e  R"  be  such  that 

Lux  +  f  (x,  ux )  <  0, 

Lu2  +  f(x,u2)>  0, 
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where  /  is  locally  Holder  continuous  in  ( x,u )  and  locally  Lipschitz  in  u  ,  and  L  is  an  elliptic 


operator  of  second  order.  Then  there  exists  a  solution  u  of  Lu  +  f(x,u )  =  0  with  ul>u>u2. 

The  Laplacian  is  a  second  order  elliptic  operator,  and  therefore  Preliminary  3-1  applies 
directly  to  equation  (1.1). 

Definition  3-2  Given  Preliminary  3-1,  we  call  u \  an  upper  solution  and  u2  a  lower  solution  to 
Lu  +  f(x,u )  =  0 . 

As  mentioned  with  the  Arzela-Ascoli  Theorem,  we  will  often  be  working  with  sequences 
of  functions.  In  many  of  our  proofs,  we  will  be  attempting  to  put  bounds  on  these  functions. 

One  useful  inequality  for  doing  so  is  GronwalTs  Inequality  or  GronwalTs  Lemma. 

Preliminary  3-2  (Theorem  1.5.7  of  [10])  (GronwalTs  Inequality)  Let  t0<t  <t{.  Let  y/{t)  and 
<f>{t)  be  continuous  functions  such  that  y/(t)  >  0  and 


<  K  +  J  y/{s)(j){s)ds 

1 0 


holds  for  t0<t<t,  and  K  a  constant.  Then 


(t  3 


(j>{t)  <  K  exp 


J  i//(s)ds 


Go  J 
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for  t0<t  <tl. 

3.2  Main  Results 

We  now  present  our  main  results  for  single  equations.  Our  first  result  is  an  extension  of 
Theorem  22  of  [26]  which  required  0  <  a  <  f3  <  1 .  We  adopt  a  similar  argument  to  obtain  results 
for  0  <  a  <  (3  <  1 . 

Theorem  3-1  Suppose  p  and  q  are  nonnegative  locally  Holder  continuous  in  M"  and 

oo  co 

[  r  max  p(x)dr  <  oo  ,  [  r  max  q(x)dr  <  oo  .  (3.1) 

J  \x\=r  J  lx|=r  v  7 

0  0 

Then  (1.1)  has  a  nonnegative  nontrivial  entire  bounded  solution  in  R"  if  0  <  a  <  f3  <  1 . 

Proof.  Since  [26]  shows  the  case  for  0  <  a  <  ft  <  1 ,  we  will  assume  0  <  a  <  [3  - 1 .  Define 

d(r)  =  m<x\{p(x),q(x)).  (3.2) 

We  will  construct  upper  and  lower  solutions  for  (1 . 1)  by  examining  the  equation 

Az  =  0{r)(za  +zp)  =  9{r)(za  +z),  r=|x|elR".  (3.3) 

Since  this  equation  is  radial,  we  can  rewrite  the  Laplacian  in  radial  fonn.  Doing  so,  equation 
(3.3)  is  equivalent  to  the  ordinary  differential  equation 
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z"(r) + 


(3.4) 


— — -z'(r)  =  0(r)(za  +z),  re[0,co). 

r 

Multiplying  each  side  of  (3.4)  by  rn_1  yields 

+  !-z'(r)  =  (r^'z'Cr))'  =  r"-10(r)(zar  +z) . 

Integrating  twice,  we  obtain 

r  t 

z(r )  =  c  +  J  f1-'1  J  1s,"_1#(.s,)(za  +  z)dsdt,  (3.5) 

0  0 

where  c  =  z(0)  >  0  is  our  central  value.  Therefore  a  solution  to  (3.3)  is  any  fixed  point  of  the 
operator  T  :  C[0,co)  — »  C[0,co)  defined  by 

r  t 

Tz{r)  =  c  +  ^tl~n^sn~ld{s){za  +z)dsdt,  0  <r  .  (3.6) 

0  0 

Note  the  integration  in  this  operator  implies  a  fixed  point  z  e  C[0,oo)  is  in  C2[0,co)  ,  making  z  a 
classical  solution  to  (3.3).  We  will  show  there  exists  such  a  fixed  point  Tz  =  z  and  that  z  is 
bounded. 

Let  z0(r)  =  c  for  all  r  >  0  ,  and  define  the  sequence 


r  t 

zk(r)  =  Tzk-\(r)  =  c  +  \  t'-n$sn-'d(s)(zl,  +  zkX)dsdt 
0  0 


for  k  - 1, 2, _  Induction  shows  this  sequence  is  increasing.  Clearly 


(3.7) 
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r  i 

z0(r)  =  c<c  +  jt1-”jsn-16>(s)(z"  +z0)dsdt  =  zx(r)  . 


Then,  assuming  zk  <  zk+1 ,  we  have 


r  i 

zk+ i(r)  =  c-\-^tl~n^sn~l6{s){zk  +zk)dsdt 
0  0 
r  t 

<  c  +  J  t1-" ^ s"~x 0{s){z%+l  +  zk+l)dsdt 
0  0 

=  zkAr)- 


Hence  {zk}  is  increasing.  We  now  show  \zk)  is  uniformly  bounded.  Integrating,  (3.7)  becomes 


r  i 

zk  (r)  =  c  + 1  tl~n  |  sn~l0{s\zak_x  +  zk_x  )dsdt 
o  o 

r  t 

<c  +  ^tl~n^sn~xO(s)(zk  +zk)dsdt 


=  c  +  - 


o  o 

1 


2  —  n 


-i  t=r 

t  r  i 

t2-n  f sn  lQ(s^z“  +  zk)ds  -  [ - t2~ntn~l0(t)(zk  +  zk)dt 

J  J  2-o 


1=0 


(3.8) 


By  L’HopitaTs  rule,  we  have 


lim 

t — >0 


l 

t1-n\sn-]d(s){z°k+zk)ds 


=  lim 

t — >0 


\sn-'d(s)(z°k  +zk)ds!tn- 
0 

=  lim[<-'6'(0(zt“+zt)/(«-2)<-!] 
=  lim[r6>(0(z“  +  z,)/(i>- 2)] 


t — >0 

=  0. 


(3.9) 


Thus  (3.8)  and  (3.9)  imply 
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zk(r)<c  + 


2  —  n 


4.1-n 


|  sn~xO(s)(zak  +zk)ds 


J  2  —  n 


t= o 


=  c  +  - 


—  r1-  J  S-Mstzi  +  Zt  )*  ■ -  f  0«>«  + )<* 

J  J  2- n 


2-n 


=  c  +  - 


n-2 


r  r 

J tO(t)(zk  +zk)dt-r2~n^sn~xO(s)(zk  +zk)ds 


<c  +  j  tO(t)(zk  +  z,  )dt. 


(3.10) 


We  will  consider  the  domain  of  z,  in  two  intervals.  Notice  that 


z[{r)  =  r'-"\sn-'d(r)(zt]+zk_,)ds  >  0  . 
0 


Thus,  we  may  choose  rk  such  that 


zk(r)<  1  for  r  e  [0,rA  ], 
z/c  (r)  >1  for  r  e  [/; ,  oo). 


(3.11) 


(3.12) 


It  is  possible  that  =  0  or  rk  =  oo  .  Indeed,  since  z[(r)  >  0  ,  we  see  from  (3.12)  and  (3.5)  that  if 
our  central  value  c>  1 ,  then  rk=  0 .  If  0  <  c  <  1 ,  then  either  rk  =  oo  with  zk  (r)  <  1  for  all  r  >  0 
or  rk  <  cc  with  zk(rk )  =  1 .  We  are  attempting  to  bound  { zk } ,  and  therefore  rk  =  oo  is  trivial. 
Therefore  if  we  consider  rk<  oo  ,  and  split  our  domain,  (3.10)  becomes 


3-6 


zk  (r)  ^  c  +  — f  ^(0(z“  +  zk  )dt  +  f  ^(0(z“  +  z*  )dt 

n-2{  n- 2J 

^  rk  \  r 

<c-\ - J  tO(t)dt  H - J  t0(t)(2zk  )dt 

n  2  0  n  ^  >k 

2  r  2  r 

<  c  ■+ - [  tO{t)dt  H - [  t8(t)zkdt 

n -2 i  n- 2 J 

0  'i 

2  1 

=  i?  -t - [  tO(t)zkdt 

n-  2J 

n 


where  R  is  a  constant  because  (3.1)  implies 


(3.13) 

<  oo. 


Finally,  Gronwall’s  inequality  (Preliminary  3-2)  yields 


zk(r)<R  + 


n  —  2 

( 


r 

J  tO(t)zkdt 


<  Rex p 


n  —  2 

v  rk 

f 


r 

J  tO{t)dt 


<  Rex p 

=  M  <  oo 


n  —  2 


J  tO(t)dt 


which  is  finite  due  to  (3.13).  We  have  shown  that  {zk}  is  a  uniformly  bounded  monotonic 
sequence.  Therefore  it  converges  pointwise  to  some  z.  Further,  since  {zk}  is  in  the  form  of 
(2.15),  Lemma  2-1  implies  that  {zk}  is  also  equicontinuous.  Pointwise  convergence  and 
equicontinuity  imply  uniform  convergence,  and  thus  z  e  C2[0,oo)  .  We  have  a  fixed  point  of 


3-7 


(3.6)  and  bounded  solution  of  (3.3).  Finally,  we  show  this  function  z  and  its  bound  M  form 
upper  and  lower  solutions  to  (1.1). 

Let  Wj  =  M  >  z  =  u2.  Clearly, 

Ai/j  =  0  <  puxa  +  qux , 

A u2  =  0(u2a  +u2)>  pu2  +qu2, 

so 

Aw,  -  (pu“  +  qux)  <  0, 

A u2  -  ( pu2  +  qu2)  >  0. 

Hence,  there  exists  a  positive  nontrivial  entire  bounded  solution  u  of  (1.1)  such  that 
M  =  ul>u>u2=z  by  the  Barrier  Method  (Preliminary  3-1).  □ 

We  will  also  consider  Theorem  2 1  in  [26]  for  0  <  a  <  /3  <  1 .  This  again  fills  the  gap  for 
when  P  - 1 . 

Theorem  3-2  Suppose  0  <  a  <  P  <  1  and  that  p(x)  =  p  (|  x  |)  e  C(R)  and  q(x)  =  q  (|  x  |)  e  C(R) 
such  that  p  and  q  are  nonnegative.  Then  equation  (1.1)  has  an  entire  large  positive  solution  if 
and  only  if 

oo  oo 

J rp(r)dr  =  oo  or  J rq(r)dr  =  oo  .  (3-14) 

o  o 
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Proof.  Again,  since  Smith  [26]  proved  this  result  for  0  <  a  <  /3  <  1 ,  we  will  consider 
0  <  a  <  P  —  1 .  An  argument  for  necessity  identical  to  that  in  Theorem  2 1  of  [26]  will  work  for 
P  —  1 .  We  are  left  to  show  sufficiency. 

Radial  solutions  of  (1.1)  satisfy  the  ordinary  differential  equation 

u"(r)  +  - — - u'(r )  =  p{\ x  | )ua  +q(\x  \ )up  .  (3.15) 

r 

As  in  Theorem  3-1,  it  follows  that  (3.15)  has  nonnegative  solutions  if  the  operator 
T :  C[0,  oo)  ->  C[0,oo)  defined  by 

r  t 

Tu{r )  =  C  +  J  t1-"  ( p(s)ua  +  q(s)u^dsdt  (3.16) 

0  0 

has  a  fixed  point  in  C[0,oo)  where  w(0)  =  c  >  0  is  our  central  value.  We  begin  by  showing  that 
(3.16)  has  a  fixed  point  in  C[0,f?]  for  arbitrary  0  <  R  <  oo  .  Similar  to  our  previous  proof,  we  let 
u0  =  c  and  define  the  sequence 

r  t 

uk+ 1  =  Tuk  (r)  =  c  +  J t1-"  jV'~'  ( p(s)uka  +  q{s)uk ) dsdt ,  (3.17) 

0  0 

k  =  1,2, _  The  same  induction  argument  from  Theorem  3-1  shows  the  sequence  {uk }  is  non¬ 

decreasing.  We  will  now  prove  that  {uk}  in  uniformly  bounded  and  equicontinuous  on 

0  <r<R. 

Define 
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0(s)  =  max  { p(s),  g(s)}  . 


(3.18) 


Integrating  as  we  did  in  (3.8)  and  (3.10),  we  have 


1 

4+1  (r)  <c  +  - — —  J  sQ{s)(uka  +  uk  )ds  ,  0  <  r  <  R 


(3.19) 


Also  similar  to  the  previous  proof,  we  have  u'k>  0 ,  so  we  may  choose  rk  such  that 


uk(r)<  1  for  r  e  [0,rA  ], 
uk (r)  >  1  for  re[rA_,f?]. 

We  are  attempting  to  obtain  a  bound,  and  thus  we  need  only  consider  rk  finite  and  r  >  rk . 
Splitting  our  integral  from  (3.19)  into  two  pieces,  we  have  by  (3.20) 


uk+l(r)  <c  +  — —  [  sQ(s)(uka  +  uk  )ds 

n~  2J0 

rk  1  r  \ 

=  c  +  [ - sQ{s)(uka  +uk)ds+  [ - s0(s)(uka  +  uk  )ds 

J  n  —  2  J  n  —  0 


n-  2 


*  2  1 

<c  +  [ - s0(s)ds+  [ - s0(s)(2uk)ds 

J  n  —  2  J  ii  -  9 


n  —  2 


c  +  J  h(s)ds  +  J  h(s)ukds 


where  h(s)  =  2s0(s)/(n-2) .  Since  {uk }  is  an  increasing  sequence  of  increasing  functions,  we 
must  have  rk+l<rt  for  all  k  —  0,1,....  Thus,  ^  for  all  k  which  yields 
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(3.21) 


'k  i 

uk+x(r )  <  c  +  J h(s)ds  +  J h(s)ukds 

0  rk 

'o  r 

<  c  +  J  h(s)ds  +  J  h(s)ukds 
0  0 
r 

=  W  +  J  h(s)ukds. 


We  now  use  induction  to  show 


*k+ 1 


r 

(r)  <  W  exp  J  h(s)ds 


u 


Clearly  we  have 


un  =  c< 


'0  /a  r 

c  +  J  h(s)ds  exp  J  h(s)ds  =  W  exp  J  h(s)ds  . 


Now,  assume  that  (3.22)  is  true  for  arbitrary  k .  Then  for  k  + 1 ,  (3.21)  implies, 


r 

uk+l  <  W  +  J h(s)uk(s)ds 
o 

<  W  +  j  h(s)W  exp  f  j  /z(t)d^ 

0  VO 

=  W  +  W  exp  j  h{t)dt 

Vo 

=  W  +  W  Qwp{^Ji(s)ds^-W 
=  W  exp  |  h(s)ds  j . 


(3.22) 
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Thus,  (3.22)  is  true  by  induction.  Finally,  since  each  uk  is  increasing,  we  arrive  at 

uk{r)<W  QxJ\^h{t)dt\^<W  QxJ\^h{t)dt^  =  M  R,  0  <r<R.  (3.23) 

That  is,  {uk}  is  uniformly  bounded  on  r  e  [0,f?] . 

Notice  that  our  sequence  is  of  the  same  fonn  as  (2.15).  Therefore  {uk}  is  also 
equicontinuous  on  r  e[0,f?]  by  Lemma  2-1.  Since  {uk}  is  a  mono  tonic,  uniformly  bounded, 
equicontinuous  sequence  of  functions  on  [0,  R],  uk  — »  u  unifonnly.  That  is  Tu  =  u  for 
re[0,i?],  and  we  have  a  fixed  point  of  (3 . 1 6)  in  C[0,  R ] . 

Next,  we  extend  this  result  to  show  that  (3.16)  has  a  fixed  point  in  C[0,oo)  .  We  do  so 
using  a  diagonal  argument  similar  to  the  argument  in  Theorem  1  of  [15].  Define  the  sequence  of 
fixed  points  { wk }  by 


Twk  =  wk  on  [0,*] ,  wk  e  C[0,&] .  (3.24) 

Restricted  on  [0,1] ,  { wk }  is  bounded  by  Mx  as  given  in  (3.23).  Using  Lemma  2-1,  we  can  also 
show  it  is  equicontinuous  on  this  interval.  Thus,  the  Arzela-Ascoli  Theorem  (Preliminary  2-1) 
implies  that  there  exists  a  subsequence,  call  it  {w\}  ,  that  converges  uniformly  on  [0,1] .  Let 

w\  — >  Vj  uniformly  on  [0,1]  as  k —>  oo  (3.25) 

Similarly,  the  sequence  {w\}  is  bounded  and  equicontinuous  on  the  interval  [0,2] .  Hence,  it 
must  contain  a  convergent  subsequence  {w2k}  that  converges  uniformly  on  [0,2].  Let 
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w\  — >  v,  uniformly  on  [0, 2]  as  k  — >  oo 


(3.26) 


Note  that  {w2k}  c=  {w*}  c:  {wt}"=2 .  This  implies  v2  =  v,  on  [0,1] .  Continuing,  we  have 
{wnk}  ci  •  •  •  c:  { vv] }  c:  !  ny  }"=n  and  a  sequence  { vy  {  such  that 

v*(r)eC[0,£]  A:  =  1,2,... 

vt(r)  =  vi(r)  for  re  [0,1] 

v*(r)  =  v2(r)  for  re  [0,2]  (3.27) 

Tt  W  =  Vi  0)  for  r  e  [0,  k- 1]. 

Thus,  we  obtain  a  sequence  { vy  ]  that  converges  to  v  on  [0,oo)  satisfying 

v(r)  =  vk(r)  if  0  <  r  <  k  .  (3.28) 

This  convergence  is  unifonn  on  bounded  sets,  implying  v  e  C[0,oo) .  We  therefore  have  our 
fixed  point  Tv  =  v  of  (3.16)  in  C[0,oo)  ,  and  equation  (1.1)  has  an  entire  radial  solution.  We 
lastly  must  show  that  this  solution  is  large. 

We  note  the  argument  of  Theorem  1  in  [15]  demonstrates  that  (3.14)  implies 

oo  t 

J  t1'"  |  .s’”-1  (p(s)  +  q(s))dsdt  =  oo .  (3.29) 

0  0 

Our  solution  satisfies 
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(3.30) 


r  t 

i/(r)  =  c  +  J tx~n^sn~x  [p(s)u*  +  q{s)u^dsdt 
0  0 

r  t 

>  c  +  min  | ua  (0), i/(0)|  J  t {~n  J sn~l  [ p(s )  +  g(s)]  dsdt 

0  0 
r  t 

=  c  +  min  |w“  (0),  i/(0)|  J  t'~n  1 [  p(s )  +  <y(.s’)] 

o  o 


Therefore,  letting  r  — >  co  ,  (3.29)  and  (3.30)  imply  lim«(r)  =  co  ,  and  we  have  a  large  solution. 

r— >oo 

This  concludes  our  proof.  □ 


The  condition  (3  =  1  remained  an  open  problem  because  comparisons  between  u"k  (r) 

and  wf  (r)  are  quite  different  depending  whether  or  not  uk(r )  <  1  or  tq(r)  >  1 .  To  resolve  this 

issue,  we  simply  split  the  domains  of  our  functions  into  two  intervals  and  carried  out  the 
comparisons  separately. 

3.3  Examples 

In  this  section  we  present  several  examples  of  equations  that  satisfy  the  restrictions  given 
in  our  analysis.  We  also  solve  several  examples  numerically.  Note  that  Theorem  3-2  only 
applies  to  radial  equations  while  Theorem  3-1  applies  to  more  general  nonradial  problems.  For 
simplicity,  all  our  examples  will  be  radial. 

Examining  Theorem  3-1,  we  must  choose  functions  p,q  that  satisfy  (3.1).  That  is, 

equation  (1.1)  must  contain  nonnegative  functions  p, q  e  C[0, oo)  that  decay  at  faster  than  Ms2 
as  s  — »  oo .  For  example, 

A«  =  (|x|+l)~V+(|x|+2)~V  or  Ay  =  2HaV  +3‘wy^  (3.31) 
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are  two  very  simple  equations  that  meet  this  requirement. 

According  to  our  analysis,  for  0  <«</?<  1 ,  entire  bounded  solutions  exist  for  equations 
of  this  form.  Choose  a  =  1/2  ,  jB  =  1 ,  and  consider 

A u  =  2 ~wu1/2  +  3~ww  .  (3.32) 

Using  a  Runge-Kutta  algorithm,  we  will  solve  this  problem  numerically.  Figure  3-1  shows  these 
solutions  for  a  variety  of  central  values  as  initial  conditions.  Indeed,  we  have  bounded  entire 
solutions. 


Figure  3-1.  Numerical  Solutions  of  (3.32)  for  Various  Central  Values 

Next,  we  consider  examples  for  Theorem  3-2.  We  now  must  have  p,g  satisfy  (3.14).  A 
few  simple  examples  would  be 
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Am  =|  x  |  ua+  |  x  |2  ,  Am  =  2|a|m"  +3wm'0  ,  Am  =  (l/(l+ 1  x  |))m“  +  (l/(l+  |  x  |))m^  , 


or 


\-2 


Am  =  (1+  |  x  |)  “  m“  +(l+ 1  x  |)  m"  . 

Again,  we  will  choose  to  examine  some  of  these  equations  numerically.  Consider 

Am  =  (l/(l+ 1  x  |))m1/2  +(l/(l+ 1  x  |))n . 


(3.33) 


(3.34) 


We  examine  a  variety  of  central  values  between  0  and  100,  and  plot  our  numerical  solutions  in 
Figure  3-2. 
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Figure  3-2.  Numerical  Solutions  of  (3.34)  for  Various  Central  Values 

We  show  results  for  0  <  |x|  <  100  and  clearly  observe  solutions  grow  quickly.  We  then  show  the 

same  solutions  in  0  <  |x|  <  5000  .  As  we  examine  larger  values  for  |x| ,  our  solutions  continue  to 

be  defined.  This  is  exactly  what  we  would  expect;  our  analysis  guarantees  these  solutions  are 
entire  and  large. 

We  will  examine  another  interesting  equation.  Consider 


Am  =  (1+  |  x  |)  "  m“  +  (1+ 1  x  |)  A 


(3.35) 
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Note  for  this  example  p(r)  =  q(r)  =  1/(1  +  r)2  satisfies  (3.14)  even  though 

lim  rp(r)  =  lim  rq(r)  =  lim  r(l  /(I  +  r)2)  =  0  .  The  integrand  decays  to  zero,  but  the  value  of  the 

r— >oo  r— >oo  r— >oo 

integral  in  (3.14)  still  approaches  infinity  as  r  — >  oo  .  Solutions  are  shown  for  a  variety  of  central 
values  in  Figure  3-3. 


Figure  3-3.  Numerical  Solutions  of  (3.35)  for  Various  Central  Values 

Our  solutions  grow  much  slower.  From  the  plot,  it  is  not  clear  that  solutions  are  large.  Flowever, 
our  analysis  guarantees  that  every  solution  goes  to  infinity  as  |  x  |— »  oo  . 
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IV.  Systems 


We  now  consider  semilinear  elliptic  systems  as  given  in  (1.2).  As  with  our  section  on 
single  equations,  we  begin  by  presenting  preliminary  material  essential  for  our  study  of  systems. 

4. 1  Preliminaries 

We  first  provide  several  definitions.  We  define  a  set  of  central  values,  what  it  means  for 
a  function  to  be  circumferentially  positive,  and  we  construct  important  functions. 

Definition  4-1  (Defined  in  Theorem  1  of  [13])  Given  system  (1.2),  we  define  the  set  of  central 
values  S  as 

S  =  {( a,b )  <e  IR+  xIR+  :  u( 0)  =  a,  v(0 )-b,  and  (u,v)  is  an  entire  solution  to  (1.2)}.  (4.1) 

Definition  4-2  (Condition  A  in  [12])  A  function  p  is  said  to  be  circumferentially  positive  (c- 
positive)  on  a  domain  Q  if  for  any  x0  e  Q.  satisfying  p(x0)  -  0,  there  exists  a  domain  Q0  such 
that  jt0efi0,  Qo  c  Q ,  and  p(x)  >  0  for  all  v  e  dCl0 . 

Definition  4-3  Given  system  ( 1 .2),  we  define 

Kl*l)  =  min[p(|x|),?(|x|)}, 

0(1  x  I)  =  max{/>(|  x  |),^(|  x  |)}.  (  ^ 

We  briefly  note  that  these  functions  have  some  useful  properties.  Namely, 

^(|  x  |)  <  p[ |  x  |)  <  0(|  x  |)  ,  ^(|  x  |)  <  <?(|  x  |)  <  <f)( |  x  |)  and  are  both  nonnegative  and  c- 
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positive  when  p,  q  are  nonnegative  and  not  identically  zero  at  infinity.  Also,  we  note  that  when 


p,q  satisfy  (2.10),  we  have 


J  ty/(t)dt  <  J  t(j){t)dt  <  J  t(p(t)  +  q(t)^dt  =  J  tp(t)dt  +  ^tq(t)dt  <  oo  . 

0  0  0  0  0 


Thus  y/,^  also  satisfy  (2.10). 


Definition  4-4  Given  system  (1.2),  we  define 


G(s)  =  min  |min /(t),  min  g(t)|,  0<s, 

H(s )  =  max  |  max  / ( t ),  max  g(t)  L  0  <  s. 

{  0 <t<s  0 <t<s  ) 


(4.3) 


These  functions  have  some  very  important  characteristics  as  well.  We  demonstrate  these 
properties  in  the  following  simple  lemma. 


Lemma  4-1  Let  f,g  e  C[0,oo)  satisfy  (2.12),  and  suppose  G  ,  shown  in  (4.3),  satisfies 
the  Keller-Osserman  condition  (2.2).  Then  G  and  H  as  given  in  (4.3)  are  continuous, 
nondecreasing,  and  satisfy 


G(s)<f(s)<H(s), 
G(s)  <  g(s)  <H(s). 


Further,  G  and  H  satisfy  (2.2)  and  (2.12). 


(4.4) 
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Proof.  It  is  clear  G  and  H  are  continuous,  nondecreasing,  and  satisfy  (2.12)  from  definition.  It 


is  also  clear  that 


G(s)<f(s)<H(s), 

G(s )  <  g(s)  <  H(s). 

Function  G  satisfies  the  Keller-Osserman  condition  (2.2)  by  hypothesis.  Then,  since 
G(s )  <  H(s ) ,  we  have 


-1/2 

dt  <  oo . 


Thus  H  satisfies  the  Keller-Osserman  condition,  and  our  proof  is  complete. 


□ 


As  we  develop  an  argument  for  the  existence  of  entire  large  solutions  for  system  ( 1 .2), 
we  will  first  show  that  solutions  exist  on  bounded  domains  for  central  values  outside  of  S  (given 
in  (4.1)).  In  order  to  prove  this,  we  will  utilize  a  version  of  Schauder  Fixed  Point  Theorem. 
Lemma  4-2  is  the  version  we  will  apply  and  is  a  very  simple  consequence  of  this  theorem. 

Preliminary  4-1  (Corollary  1 1 .2  of  [9])  (Corollary  of  Schauder’s  Fixed  Point  Theorem)  Let  G 
be  a  closed  convex  set  in  a  Banach  space  B  ,  and  let  T  be  a  continuous  transformation  of  G  into 
itself  such  that  the  image  TG  is  precompact.  Then  T  has  a  fixed  point. 


Definition  4-5  A  linear  transformation  between  two  Banach  spaces  is  called  compact  (sometimes 
called  completely  continuous)  if  the  images  of  bounded  sets  are  precompact. 
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Preliminary  4-2  (Theorem  5.24.2  of  [19])  Any  compact  linear  transformation  between  Banach 
spaces  is  continuous. 

Lemma  4-2  Let  G  be  a  closed,  convex,  and  bounded  set  in  a  Banach  space  B  ,  and  let  T  be  a 
compact  transformation  that  maps  G  into  itself.  Then  T  has  a  fixed  point. 

Proof.  This  follows  directly  from  Preliminary  4-1,  Definition  4-5,  and  Preliminary  4-2.  □ 

After  using  Lemma  4-2  to  establish  solutions  on  bounded  domains,  we  will  show  that  any 
bounded  solution  on  such  a  domain  may  be  extended  to  include  a  larger  domain.  We  accomplish 
this  extension  with  Cartheodory’s  Theorem.  We  now  provide  this  theorem  and  related 
definitions. 

Definition  4-6  Consider  the  equation 


x'  =  f(t,x)  (4.6) 

with  initial  condition  x(z)  =  %  .  We  say  that  ^  is  a  solution  to  (4.6)  in  the  extended  sense  on 
some  interval  /  if  <p'(t )  =  f{t,<f>{t ))  for  all  ie/  except  on  a  set  of  Lebesgue -measure  zero  and 
¥j)  =  £  • 

Preliminary  4-3  (Theorem  1.1  of  Chapter  2  of  [7])  (Caratheodory’s  Theorem)  Consider  the 
equation  given  in  (4.6).  Let  a,b  be  constants,  and  define  the  rectangle 
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R  =  j(fx)  :  \t-r\  <  a,  |x-^|<b}  where  (r,£)  is  a  fixed  point  in  the  ( t,x )  plane.  Let  /  be 
defined  on  R  .  Suppose  /  is  measurable  in  t  for  each  fixed  x ,  and  continuous  in  x  for  each 
fixed  t .  If  there  exists  a  Lebesgue-integrable  function  m  on  the  interval  \t  -  r|  <  a  such  that 

\f{t,x)\<m{t),  ( t,x)eR , 

then  there  exists  a  solution  (j)  of  (4.6)  in  the  extended  sense  on  some  interval  |  t-r  \<J3,(J3>0) 
satisfying  (/){ r)  =  % . 

Note  that  if  /  is  continuous  on  R  ,  then  ^  is  a  solution  in  the  traditional  sense.  That  is 
</>'{()  is  continuous  and  (/. )  =  / (t,c/>(t  j)  for  all  t  e  (t  -  J3,t  +  J3) .  We  are  now  ready  to  present 
our  main  results  for  systems. 

4.2  Main  Results 

The  theorems  and  lemmas  presented  in  this  section  build  off  one  another.  We  begin  by 
proving  existence  of  entire  solutions  and  work  toward  showing  the  existence  of  entire  large 
solutions. 

Theorem  4-1  Assume  p,q  e  C(R")  are  nonnegative,  not  identically  zero  at  infinity,  and  satisfy 
(2.10).  Let  /,ge  C[0,oo)  satisfy  (2.12),  and  let  G  ,  given  in  (4.3),  satisfy  the  Keller-Ossennan 
condition  (2.2).  Then  system  (1.2)  has  infinitely  many  entire  nonnegative  solutions. 

Proof.  Radial  solutions  of  system  (1.2)  are  solutions  to  the  ordinary  differential  system 
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u"(r)  +  ^-^-u'(r)  =  p(r)f  (v(r)),  r  >  0, 

r 

v"(r)  +  — — -v'(r)  =  q(r)g(u(r)),  r  >  0. 

r 


It  follows  that  solutions  to  ( 1 .2)  are  fixed  points  of  the  operator 
T(u,v) :  C[0, oo) x  C[0, oo)  ->  C[ 0, oo) x  C[0, oo)  defined  by 


T(u(r),v(r))  =  - 


A  l 

a  + 1  tx~n  J  sn~lp(s)  f  ( v(s))dsdt ,  r  >  0, 
0  0 
r  t 

b  +  J  f"  J  snlq(s)g(u(s))dsdt,  r>  0 


where  u(  0)  -a>  0  and  v(0)  =  b>  0  are  the  central  values  for  the  system.  We  will  begin  by 
establishing  a  fixed  point  of  (4.8)  in  C[0,f?]xC[0,f?]  for  arbitrary  R  >  0  . 

Define  u0  (r)  -  a  and  v0(r)  =  b  for  all  r  >  0 ,  and  consider  the  sequence 


(uk,vk)  =  T(uk_  i,v*_1)  =  - 


a  +  J  t]  s'1  1  p(s)  f  (vk_ j  ( s))dsdt ,  r  >  0, 

0  0 
r  t 

b  +  J  f J s',^q(s)g(uk_l  ( s))dsdt ,  r  >  0 


for  A:  =  1, 2 ... .  Consider  the  single  equation  u  given  by 


A u  =(/>{\x  | )//(«),  x  i 


m  — >  oo  as  x  — >  oo 


(4.10) 
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where  (f>  and  H  are  defined  in  (4.2)  and  (4.3)  respectively.  Recall  that  (j)  is  nonnegative,  c- 
positive,  and  satisfies  (2.10).  By  Lemma  4-1,  H  satisfies  the  Keller-Osserman  condition  (2.2). 
These  conditions  guarantee  such  a  17  exists  by  Theorem  2  of  [12],  Further,  this  radial  solution 
satisfies 


r  t 

u(r )  =  17(0)  +  jY“"  J  sn~l 0(s)H (u(s))dsdt ,  r>  0 
0  0 


where  u  (0)  >  0 .  Note  this  solution  is  non-decreasing.  Indeed 


(4.11) 


—  u(r )  =  r1  ' 
dr 


1 


sn  lt/>(s)H(u(s))ds  >  0 . 


(4.12) 


We  shall  choose  our  central  values  a  and  b  such  that  0  <  a  <  u  (0)  and  0  <  b  <  u  (0) .  Next,  we 
use  induction  to  show  each  individual  sequence  in  (4.9)  is  uniformly  bounded  on  bounded  sets. 
Specifically,  we  will  show  ma x{uk(r),vk(r)}<u(R)  for  all  k  and  all  re[0,R]. 

Since  H  and  u  are  non-decreasing  functions  and  a,b<  «(()) ,  we  have 


r  t 

ux{r)-  a  + 1 tx~n | / (v0 (s))dsdt 
0  0 

r  t 

<  17(0)  +  J  tl-n  J  snl</>(s)H(b)dsdt 

0  0 
r  t 

<  17(0)  +  J  tl^snl</>(s)H(u(s))dsdt 

0  0 


=  U  (r), 

and  similarly  v, (r)  <  u(r).  Now,  assume  uk ,  vk  <  17  .  This  implies 
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r  t 

uk+ 1  (r)  =  a  +  J  t^n  J  sn~xp(s)f(vk  (, s))dsdt 
0  0 

r  t 

<  u  (0)  +  J  tx ""  J  sn~x(j)(s)H{yk  ( s))dsdt 

0  0 
r  t 

<  w(0)  +  JY‘"  J sn  X<f>(s)H(u(s))dsdt 

0  0 

=  u(r). 


We  similarly  obtain  vk+1(r)  <  u(r),  so  uk,vk  <u(r)  by  induction.  Again,  since  u  is  non¬ 
decreasing,  we  have  the  bound 


uk(,r),vk{r)<u{r) 

<u(R) .  (4.13) 

=  M 

for  all  k  >  1  and  all  r  e  [0,  f?] .  We  know  M  <  oo  because  u  is  entire. 

Then,  since  each  individual  sequence  {uk},{vk}  is  of  the  form  (2.15),  Lemma  2-1 
guarantees  that  {uk},{vk}  are  each  equicontinuous  on  [0,f?] .  We  have  two  sequences,  {uk }  and 
{ vk } ,  that  are  uniformly  bounded  and  equicontinuous  on  [0,  R\ .  By  the  Arzela-Ascoli  Theorem 
(Preliminary  2-1)  there  exists  a  u  e  C[0,  R\  and  a  subsequence  {uk  }  such  that 


uk  — >  u  unifonnly  on  [0,  R\ . 


If  we  consider  the  corresponding  subsequence  {vk  }  c  [v/(  J ,  the  Arzela-Ascoli  Theorem  again 
implies  there  exists  ave  C[0, R\  and  a  subsequence  of  {vA  } ,  call  it  {vk  }  ,  such  that 

j  Ji 


vk  —>  v  uniformly  on  [0,f?]. 

Ji 
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Therefore  we  have 


(uk  ,vk  )  — >  (w,v)  unifonnly  on  [0,f?]x[0,f?] . 

Ji  Ji 

Hence,  (u,v)  is  a  fixed  point  of  (4.8)  in  C[0,f?]x  C[0,f?] .  Next,  we  extend  this  result  to  show  T 
has  a  fixed  point  in  C[0,  qo)  x  C[0,  qo)  .  We  proceed  with  a  diagonal  argument  similar  to  that  in 
Theorem  3-2.  However,  in  this  case,  we  are  considering  fixed  points  and  convergence  in 
C[0,  qo)  x  C[0,  qo)  rather  than  C[0,oo)  .  The  idea  is  the  same. 

Let  {(wk,zk)}  be  a  sequence  of  fixed  points  defined  by 

T(wk,zk)  =  (wk,zk)  on  [0,&],  (wt,zt)eC[0,/c]xC[0,/c],  (4.14) 

for  k  =  1,2,... .  As  earlier,  we  may  show  that  both  { wk }  and  { zk }  are  bounded  and 
equicontinuous  on  [0, 1] .  Thus  by  applying  the  Arzela-Ascoli  Theorem  to  each  sequence 
separately,  we  can  derive  {(wk,zk)}  contains  a  convergent  subsequence,  {(w\,z\)} ,  that 
converges  uniformly  on  [0, 1]  x  [0, 1] .  Let 

(w\,z\)  — »  (Wj,Vj)  unifonnly  on  [0, 1] x [0, 1]  as  k  — >  oo  . 

Likewise,  the  subsequences  { vv,' }  and  {z\}  are  each  bounded  and  equicontinuous  on  [0,2]  so 
there  exists  a  subsequence  {( w2k,z2k )}  of  {(w\,z\)}  such  that 

(w2k,z2k)  — >  (m2,v2)  unifonnly  on  [0,2]x[0,2]  as  k —> qo  . 
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Notice  that  {(w2k ,z\)}  c  {(w\,z\)}  c  {(w{,z*)}”  2  so  (n2,  v2)  =  (w,,  v,)  on  [0,l]x[0,l]. 

Continuing,  we  obtain  a  sequence  {(14,^)}  such  that 

(n,,v,)eC[0,£]xC[0,£] 

(Mi(r),vA.(r))  =  (w1(r),v1(r)) 

(M*0*),vt(r))  =  (M2(r),v2(r)) 

(«*  W.  fit  W)  =  («*_i  (r),  vw  (r)) 

Thus  (uk,vk)  converges  to  (u,  v)  which  satisfies 

(u(r), v(r))  =  (uk (r), vt (r))  if  0 <r<k  (4.15) 

The  convergence  is  unifonn  on  bounded  sets,  and  thus  (u(r),v(r))  e  C[0,  qo)  x  C[0,  qo)  is  a  fixed 
point  of  (4.8)  and  an  entire  solution  to  (1.2).  We  chose  our  central  values  0  <  a  <  w(0)  and 
0<b<w(0)  arbitrarily  where  u  is  defined  in  (4.10).  Therefore  [0,i7(0)]  x[0,i7(0)]  is  a  subset  of 
our  set  of  central  values  S  given  in  (4.1).  We  conclude  (1.2)  has  an  infinitely  many  entire 
solutions.  □ 

Note  we  did  not  use  the  function  G  from  in  (4.3)  directly  in  this  result.  Rather,  since  G 
satisfied  the  Keller-Osserman  condition  (2.2),  we  were  guaranteed  H  satisfied  this  condition 
according  to  Lemma  4-1.  We  then  used  H  in  our  argument.  The  function  G  need  not  satisfy 
(2.2)  to  show  existence  of  entire  solutions.  In  fact,  we  only  need  /  or  g  to  satisfy  the  Keller- 
Osserman  condition,  which  we  prove  in  the  following  corollary. 


k  =  1,2,... 

for  r  e  [0, 1] 
for  r  e  [0,2] 

for  r  e[0,£-l] 
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Corollary  4-1  Assume  p,q  e  C(M")  ,  are  nonnegative,  not  identically  zero  at  infinity,  and  satisfy 
(2.10).  Let  f,geC[ 0,oo)  satisfy  (2.12).  Also,  suppose  /  or  g  satisfies  the  Keller-0 ssennan 
condition  (2.2).  Then  system  (1.2)  has  infinitely  many  entire  nonnegative  solutions. 


Proof.  Without  loss  of  generality,  suppose  /  satisfies  (2.2).  Then  H  as  defined  in  (4.3)  still 
satisfies  the  Keller-Ossennan  condition  because 


-1/2 


J  J  H{s)ds  dt  <|  j  / ( s)ds 


-1-1/2 


dt  <  co. 


In  addition,  H  clearly  must  still  satisfy  (2.12).  Our  argument  then  proceeds  exactly  like  the 
proof  of  Theorem  4- 1 .  □ 


While  we  may  have  the  existence  of  entire  solutions  to  system  (1.2)  when  G  does  not 
satisfy  (2.2),  it  is  necessary  for  G  to  satisfy  the  Keller-Ossennan  condition  as  we  characterize 
our  set  of  central  values.  We  examine  this  set  next. 


Theorem  4-2  Given  the  hypotheses  in  Theorem  4-1,  the  set  of  central  values  S ,  given  in  (4.1),  is 
closed  and  bounded. 


Proof.  To  show  S  is  bounded,  we  proceed  in  a  manner  comparable  to  Cirstea  and  Radulescu  [6]. 
For  contradiction,  suppose  S  is  unbounded.  Consider  the  equation 
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iel", 


(4.16) 


A77  =  ^(|x|)G(77/2), 

77  — >  oo  as  |  x  |— >  oo 

where  yr  and  G  are  given  in  (4.2)  and  (4.3).  Recall  y/  is  nonnegative,  c-positive,  and  satisfies 
(2.10).  Also,  G  satisfies  (2.2),  so  an  entire  large  solution  rj  exists  for  equation  (4.16)  by 
Theorem  2  in  [12].  Since  we  are  assuming  S  is  unbounded,  we  can  find  central  values  a, b  e  S 
such  that  a  +  b>  ?](0) .  Notice  that 


f(y(r))>G(y(r))>G 

g(u(r))>G(u(r))>G 


^u(r)  +  v(r)^ 

V  2  J 

'  u(r)  +  v(r)^ 
{  2  y 


if  v(r)  >  u(r), 
if  u(r)  >  v(r). 


Then  using  (4.17),  we  have 


(4.17) 


A  (u  +  v)  =  p(r  )/(v)  +  q(r)g(u ) 
^V(r)(f(v)  +  g(u)) 

>V,(r)G(“(r)  +  V(''}'' 


(4.18) 


Consider  some  closed  finite  ball  B(0,R) .  We  will  use  the  maximum  principle  to  show  u  +  v<rj 
in  this  ball. 

Since  rj  and  u  +  v  are  radial,  consider  these  equations  in  their  related  ordinary 
differential  equation  form.  Define  h(r)  =  (1  +  r2)  1/2 ,  and  suppose  for  contradiction  that 
u+v>  rj  at  some  point  in  [0,7?] .  Let  6  >  0  be  small  enough  such  that 


max  [ (u  +  v)(r)  -  r/(r )  -  sh{r)  1  >  0  . 

re[0,/?]L  J 
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Let  r0  e  [0,  R]  be  the  point  where  this  maximum  occurs  so  (u  +  v)(r0)  -  77  (r0)  -  sh(r0)  >  0 . 
Therefore  at  r0 ,  by  (4. 1 8)  and  our  assumption  that  w(r0)  +  v(r0)  >  r/(r0) ,  we  have 


0  >  A ((u  +  v)(r0)  -  7)(r0)  - sh(r0)) 


>  i//(r0)G 


u(r0)  +  v(r0) 


~W( r0)G(rj/2)-A£h(r0 ) 


=  V(r0) 

>  - sAh(r0 ) 

>0, 


u(r0)  +  v(r0) 


-G(tj/  2) 


-A  eh(r0) 


(4.19) 


where  the  last  inequality  holds  because  Ah(r)  <  0  for  all  r  e  [0,f?] .  This  may  be  seen  by  direct 
calculation,  and  is  also  a  specific  fonn  of  Lemma  2.0.18  of  Proano  [23].  We  have  a 
contradiction,  and  thus 


u  +  v  <  r]  in  B(0,R) .  (4.20) 

However,  this  now  contradicts  our  original  assumption  that  a,b  were  chosen  such  that 

i/(0)  +  v(0)  =  a  +  b>r/(  0) . 

Therefore,  we  conclude  that  S  is  bounded. 

Now  we  prove  S  is  closed  by  showing  S  contains  its  boundary.  Let  (a0,b0)  e  dS  .  Then 
there  exists  some  ball  centered  at  ( a0,b0 )  with  radius  1  Ik  such  that  B((a0,b0),  1/ k) f)  S  ^  0 . 

For  each  k  >  1 ,  we  denote  the  arbitrary  point  (ak0,b^)  e  ST) B((a0,b0),  1/ k ) .  Note 
{(ak0  ,bk0 )}  — >  (a0,b0)  as  k  —>  co  .  Then,  we  define  the  sequence 
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(4.21) 


ao  +  J ^  s ”  1  p(s)  f  (vk  ( s))dsdt ,  r  >  0, 

(uk(r),vk(r))  =  <  °  ® 

bo  +  J  b  " J  s"  l(l(s)g(.uk  ( s))dsdt ,  r>  0 

o  o 

where  each  (uk ,  vk )  is  an  entire  solution  to  (1.2).  These  solutions  exist  since  each  central  value 
(ak0,bl)  e  B((a0,b0),  1/ k) f|  S  a  S  .  We  now  show  {(uk,vk)}  has  a  convergent  subsequence  on 
C[0,  oo)  x  C[0,  oo)  .  Similar  to  our  argument  for  Theorem  4- 1 ,  we  first  demonstrate  that  the 
sequence  has  a  convergent  subsequence  on  C[0,f?]x  C[0,f?]  for  arbitrary  R  ,  and  then  we  extend 
to  C[0,  oo)  x  C[0,  oo)  .  We  must  make  several  minor  adjustments  to  achieve  our  bounds  for  this 
new  sequence,  but  the  idea  is  comparable. 

For  each  k  =  1,2...,  (4.20)  gives  us  uk  +vi-  -  77  for  r  e  [0,  f?]  where  77  is  defined  in 
(4.16).  Since  Jj'(r)  >  0  ,  we  have 

uk(r)  +  vk(r)<rj(R)<co  for0<7 ~<R  (4.22) 

where  r/(R)  is  finite  because  77  is  entire.  Thus  {uk},{vk}  are  each  uniformly  bounded  on  [0,i?]. 
Lemma  2-1  implies  {uk},{vk}  are  equicontinuous  on  [0,  R] .  As  in  the  previous  proof,  we  may 
use  the  Arzela-Ascoli  Theorem  to  show  a  convergent  subsequence  of  (4.21)  exists  on 
C[0,f?]x  C[0,f?] .  This  gives  us  a  solution  to  system  (1.2)  in  B(0,  R)  a  M"  with  central  values 
( a0,b0 )  e  GS  .  Since  R  is  arbitrary,  we  can  use  the  same  diagonal  argument  from  Theorem  4-1  to 
show  (4.2 1)  has  a  convergent  subsequence  on  C[0, 00)  x  C[0, 00)  .  Therefore  we  have  the  solution 
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r  t 

u(r)  =  a0  +  ^tl~n  ^  sn~l  p(s)f  (v(s))dsdt9  r  >  0, 

0  0 

< 

r  t 

v( r)  =  b0  +  J  tl~n  | sn  lq(s)g(u(s))dsdt,  r  >  0 
0  0 

where  ( a0,b0 )  e  dS  .  Hence  ( a0,b0 )  e  S  implying  S  is  closed.  This  completes  our  proof.  □ 

Thus  far,  we  have  established  the  existence  of  entire  solutions  to  system  ( 1 .2)  and 
characterized  the  set  of  central  values  S  as  closed  and  bounded.  For  monotonic  and  convex 
functions  /  and  g ,  Peng  and  Song  [22]  further  characterized  S  as  convex.  We  cannot  show 
this  for  our  non-monotone  problem.  In  fact,  we  will  provide  an  example  later  for  which  S 
(found  numerically)  appears  non-convex.  Still,  Theorem  4-1  and  Theorem  4-2  provide  a  limited 
geometric  description  of  S  which  we  present  as  the  following  corollary. 

Corollary  4-2  Given  the  hypotheses  in  Theorem  4-1,  our  set  of  central  values  S  given  in  (4. 1), 
satisfies  c  5  c  f  where  7]  =  {(a,b) :  0  <a<  w(0),  0  <b<  w(0)}  and 
T2  =  {( a,b ) :  a  +b  <  77(0)}  .  The  functions  u  and  77  are  given  in  (4.10),  and  (4.16)  respectively. 
See  Figure  4-1. 
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v(0)  =  6 


Figure  4-1.  Entire  Solution  Existence  Region  for  Central  Values  of  (1.2) 


Proof.  In  Theorem  4-1,  we  proved  the  existence  of  entire  solutions  to  (1.2)  by  constructing 
Tx  c=  S  .  In  Theorem  4-2,  and  we  demonstrated  S  to  be  bounded  by  constructing  T2ziS.  □ 


Note  that  we  drew  S  in  Figure  4-1  as  connected,  but  we  have  not  shown  this  to  be  the 
case.  Indeed,  we  know  much  less  about  our  set  of  central  values  as  compared  to  monotonic 
problems.  However,  the  fact  that  S  is  closed  and  bounded  is  enough  to  show  the  existence  of 
entire  large  solutions.  Before  we  proceed,  we  will  need  to  prove  several  minor  lemmas. 
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Lemma  4-3  Under  the  hypotheses  of  Theorem  4-1  and  for  any  ( c,d )  e  K+  xffif  such  that 
(c,  d)  £  S  and  c  ^  0  ^  d  ,  system  ( 1 .2)  has  a  solution  in  some  ball  B( 0,  p ) ,  where  0  <  p  <  oo  . 

Proof.  Let  ( c,d )  be  defined  as  above.  We  wish  to  find  a  radial  solution  to  (1.2)  in  B(0,p), 
p  >  0  .  This  solution  is  a  fixed  point  of  the  operator  T  :  C[0,  p\  x  C[0,  p\  — >  C[0,  p\  x  C[0,  p\ 
defined  by 


T(u(r),v(r ))  =  (w(r),v(r))  = 


t 

c  +  j 

V""j 

C 

l  0 

r  t 

d  + 

[<"■/ 

__ 

0  0 

(4.23) 


i  <  r  <  p. 


We  will  show  for  p  sufficiently  small,  a  fixed  point,  and  thereby  a  solution  exists.  We 
accomplish  this  using  the  version  of  Schauder’s  Fixed  Point  Theorem  given  in  Lemma  4-2. 

First,  we  establish  a  subset  X  cz  C[0,  p\  x  C[0,  p\  that  satisfies  the  necessary  hypotheses  of  this 

lemma.  Note  that  for  p  >  0  ,  C[0,p]xC[0,p]  is  a  Banach  space  with  norm 

||(m,v)||op  =max(||M|oo,|v|M  where  ||w||^  =  sup  \u{r)\ .  Define  the  subset  IcC[0,/)]xC[0,/i]  by 

0<r<p 


X  =  |(w,v)  e  C[0,p]xC[0,/)] :  |(M,v)-(c,J)||aj  <  min{c,r/}| , 

where  c,  d  are  our  given  central  values.  Since  this  is  a  closed  ball  in  a  Banach  space,  X  is 
closed,  bounded,  and  convex 

We  turn  our  attention  to  the  operator  T  given  in  (4.23).  We  show  T  is  a  compact 
operator.  Let  Y  be  an  arbitrary  bounded  set  given  by 
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Y  =  |(w,v)  e  C[ 0,  >c]  x  C[0,  p] :  ||(w,  v)|oo  <  m} 


for  some  M  >0  .  Then,  suppose  (uk,vk)  is  an  arbitrary  sequence  in  Y,  and  consider  the  image 
of  this  sequence  T  (uk  ,vk)  =  (uk  ,vk) .  Clearly 


r  t 

uk  (r)  =  c  +  J  t1"’  J  sn~lp(s)f(vk  (, s))dsdt 
0  0 
r  t 

<  c  +  jY~"  J s'l  l(j){s)H{vk{s))dsdt 

0  0 
p  I 

<  c  +  jY~"  J  s"  lt/>(s)H (M)dsdt 

0  0 


where  (j)  and  H  are  given  in  (4.2)  and  (4.3),  and  M0  is  constant.  Similarly  vk(r)  <M0 .  Thus, 
uk  and  vk  are  uniformly  bounded.  Further,  since  uk(r)<M  and  vk(r)<M  for  all  k  and  all 
0  <r  <  p ,  Lemma  2-1  implies  that  uk  and  vk  are  also  equicontinuous  on  0  <  r  <  p  .  Finally, 
using  the  Arzela-Ascoli  Theorem  as  before,  there  exists  a  uniformly  convergent  subsequence 
(uk  ,vk  )  — >  (u,v)  where  (u,v)  e  C[0,p]x  C[0,p] .  This  implies  that  the  image  of  an  arbitrary 

bounded  set  under  T ,  is  (sequentially)  compact.  That  is,  T  is  a  compact  operator. 

Lastly,  we  must  show  that  T  maps  elements  of  X  back  into  X  .  Again,  take  any 
(u,v)  e  X  .  Since  X  is  bounded,  suppose  ||(m,v)||oo  <Q .  By  integrating  as  in  (3.8)  and  (3.10), 
we  have  the  estimate 
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0  0 


r  i  r  i 

J  tx~n  ^  sn~x  p{s) / ( v(s))dsdt  <  J  tx~n  J  sn~x  (j)(s)H  {y(s))dsdt 

0  0 

r  t 

\  //(v(r))J  tx~n^sn~l(/){s)dsdt 
o  o 

r 

;  H(Q )J  t<p(t)dt 
0 

=  HiQ)p\  </>(t)dt 


< 


< 


where  </>  and  H  are  given  in  (4.2)  and  (4.3).  Similarly,  we  may  show 


r  i  p 

J  tX  n^ sn~lq(s)g(u(s))dsdt  <  H(Q)p | </>(t)dt . 

0  0  0 


Since  <j)(r )  is  well  defined  for  r  e  [0,oo)  and  c,d  >  0  by  hypothesis,  we  can  choose  p>  0  small 

p 

enough  so  that  H(Q)p J  <f>(t)dt  <  min  {c,d}  .  Doing  so,  for  (w,v)  el,we  have 


T(u(r),v(r ))  =  (j?(r),v(r))  where 


c  <  u(r) 

r  t 

=  c  + 1 J  snXp(s)f{v{s))dsdt 

0  0 

p 

<c  +  H(Q)p J  (f)(t)dt 
o 

<  c  +  mm{c,d). 


and  similarly  d  <  v(r )  <  d  +  min{c,  d}  for  all  0  <r  <  p  .  Thus  (u,  v)  e  X  . 

We  have  shown  X  to  be  closed,  convex,  and  bounded.  Also,  we  have  proven  the 
operator  T  in  (4.23)  is  compact,  and  we  have  shown  for  any  x  e  X  ,  Tx  e  X  .  It  follows  from 
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Lemma  4-2  that  (4.23)  has  a  fixed  point.  Hence  system  (1.2)  has  a  solution  in  B(0,p) .  Further, 


we  chose  (c,d)&S  implying  p<  oo.  □ 

Now  that  we  have  established  the  existence  of  solutions  for  any  pair  of  positive  central 
values,  we  will  show  for  any  such  (c,  d)  £  S  ,  a  large  solution  exists  on  a  finite  domain. 

Lemma  4-4  Given  the  hypotheses  of  Theorem  4-1,  let  ( u,v )  be  a  solution  to  (1.2)  with  central 
values  (c,d)<£S  ,  c  ^  0  ^  d  ,  and  define  the  set 

Rso i  =  {r  >  0 :  there  exists  a  solution  of  (1.2)  in  5(0, r)  such  that  (i<(0),v(0))  =  ( c,d )}  .  (4.24) 

Let  Rc  d  be  given  as 


Rc  d  =  SUP  R,ol  ■ 


(4.25) 


Then  lim  u(r)  =  co  =  lim  v(r)  . 

r  ~^Rc,d  r  ~^Rc,d 

Proof.  Take  (c,d)  £  S  ,  c  ^  0  ^  d  .  By  Lemma  4-3,  Rsol  ^  0  ,  and  since  (c,d)iS  ,  Rc  d  <  oo  . 

Let  (w,v)  be  a  solution  of  (1.2)  in  B(0,Rcd )  with  central  values  ( c,d ) .  Since  u'  >  0  and  v'  >  0, 
lim  u{r)  and  lim  v(/")  exist  (possibly  infinity).  For  contradiction,  suppose  lim  u(r)  =  A  <  co  . 

r^RC,d  r^RC,d  r^Rc,d 

This  implies 
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r  i 

lim  v(r)  =  lim  d  +  [  tl~n  [  sn  l q(s)g(u(s))dsdt 

r—>R„  H  r—>R„  ,  J  J 

0  0 

Km  t 

<d+  J  tl~n ^ sn~1q(s)H(u(sJ)dsdt 


o  o 
r,.„ 


< 


d  +  H(A )  J  tn\sn-'q(s)dsdt 


0  0 


=  B  <  CO. 


Thus  u(Rcd)  =  A  and  v(Rcd)  =  B  are  well  defined.  We  now  consider  system  ( 1 .2)  over  the 
interval  Rcd<r<Rcd+s,  where  s  >  0 .  Again,  our  equations  are  radial,  so  we  may  integrate 
system  (1.2)  to  obtain 


r  t 

u(r)  =  A  +  J  f 1_" | s"~lp(s)f(v(s))dsdt,  Rcd  <r<  Rcd  +  s, 

Kd  o 

r'  ,  (4-26) 

v(r)  =  B  +  J  tx~n | sn  lq(s)g(u(s))dsdt,  Rcd  <r<  Rcd  +  s. 

Km  0 


This  is  equivalent  to  the  problem 


u\r)  =  r1_"|  sn~l  p(s)f(v(s))ds,  u(Rcd )  =  A, 

0 

r  t 

v(r)  =  B+  J  tl~"  j  sn"1q(s)g(u(s))dsdt. 

Km  0 


Rcid^r^Rcj+e> 

Rc,d  -  r-  Rc,d  +  £- 


(4.27) 


Substituting  v(r)  into  our  equation  for  u\r) ,  we  get  the  initial  value  problem 
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r 

u\r )  =  r  n^sn  lp{s)f{y(s))ds 
0 

r  st 

=  rl~"  J  sn  lp(s)f  B+  J  tx~n  |  ynXq(y)g(u(y))dydt 

RcJ  L  0 

=  F(r,u(r )) 

where  Rc  ,d^r^Rc,d+£  and  u(Rcd  )  =  A  .  If  we  show  this  initial  value  problem  has  a  solution 
on  some  interval  Rcd<r<Rcd+f3  (0  <  J3  <e),  then  (4.27)  and  (4.26)  will  have  a  solution  as 
well.  This  would  imply  that  system  ( 1 .2)  has  a  solution  on  0  <r  <Rcd  +  J3 ,  contradicting  the 
definition  of  Rc  d  in  (4.25).  We  use  Caratheodory’s  Theorem  shown  in  Preliminary  4-3. 

Define  the  rectangle 

R  =  j(r,w)  :  |r-7?C(i|  <  e,  |w-^4|<v4j  (4.29) 

where  £  >  0 ,  A  =  lim  u(r),  and  R  ,  is  given  in  (4.25).  Treating  r  and  u  as  two  independent 

r^RcA 

variables,  we  see 


i  ■>  i 

F(r,u)  =  rx~n  |  sn^p(s)f  B  +  J  tx~n  J  yn^q(y)g(u)dydt  ds 


Rc4  0 


r  hi 

=  r'  "  J  s"  lp(s)f  B  +  g{u)  J  tl"jynlq(y)dydt  ds. 


K.d  o 


(4.30) 


Recall  /,  g,  p,  q  are  all  defined  on  [0,  qo)  ,  and  consider  F(r,  u )  on  [Rc  d  -  /?,  Rc  d  +  p\  x  [0, 2  A] . 
Let  N  =  max  g(u) ,  and  choose  p  small  enough  so  that 

0<u<2A 
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RC,d~P  t  RcJ  t 

g(u)  J  t1~njy’~lg(y)(fydt  =  -g(u)  J  t'~n  ^  yn' q(y)dydt 

Rc,d  0  RC,J-P  0 

ReJ  t 

>-iV  J  t'~n^yn~'q(y)dydt 

RcJ-P  0 

>-B. 


This  ensures  the  argument  of  /  remains  positive,  and  hence  F(r,u )  is  defined  on  all 
\-Rc,d~P’Rc.d  +  p]x[0,2A] .  Then,  choosing  s  <  p ,  (4.29)  implies 

R  c=  [f?c  d  -  p,  Rc  d  +  p\  x  [0, 2^4] .  Also,  since  R  is  closed  and  bounded,  ma x  F(r,u)  =  M  is  well 

1  ( r,u)eR 

defined,  and  we  have 


\F(r,u)\  =  F(r,u)<M  for  all  (r,u)  e  R .  (4.31) 

Applying  Cartheodory’s  Theorem  (Preliminary  4-3),  (4.28)  has  a  solution  in 
Rcd-j3<r<Rcd+j3  (0  <  J3<e)  with  u(Rcd)  =  A  .  We  have  our  contradiction,  and  therefore 

lim  u{r)  =  co  .  Similarly,  we  can  obtain  lim  v(r)  =  oo ,  and  our  proof  is  complete.  □ 

r^>Rc,J  r->RCJt 


We  finally  have  the  necessary  tools  to  present  our  results  for  existence  of  entire  large 
solutions  for  systems.  Our  general  argument  is  taken  from  Lair  and  Wood  [13]. 

Theorem  4-3  Given  the  set  of  central  values  S  in  (4.1),  let 

E  =  {(a,b)  e  R+  xM+  :  (a,b)  e  dS,  0,  6^0}.  (4.32) 
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Under  the  hypotheses  of  Theorem  4-1,  any  positive  radial  solution  of  system  (1.2)  with  central 


values  ( a ,  b)sE  is  an  entire  large  solution. 


Proof.  Let  ( a,b )  e  E .  We  have  shown  in  Theorem  4-2  that  SS  cz  S  ,  so  define  (u,  v)  to  be  an 
entire  solution  with  central  values  ( a,b ) .  Since  ( a,b )  edS  ,  we  may  choose  a  point  ( ak,bk)<£  S  , 
ak  >  0 ,bk  >  0  ,  where  ( ak,bk )  e  B((a,b),\/ k) .  From  Lemma  4-3,  there  exists  a  solution  to  (1.2) 
in  the  finite  ball  5(0,  Rn  b  )  with  central  value  ( ak,bk ) .  For  simplicity,  we  will  abbreviate 
R  =  R  .  We  write  each  solution  as 

ak  A  k 


r  t 

uk  (r)  =  ak  +  J  t1"’  |  snlp(s)f{yk  ( s))dsdt , 
0  0 
r  t 

vk(r)  =  bk  +  J  f1-"  JV' lq(s)g(uk (s))dsdt, 
0  0 


re[0,Rk), 

re[0,Rk). 


From  Lemma  4-4,  we  have  for  each  Rk , 


lim  uk  (r)  =  lim  v,  (r)  =  oo  . 

r^Rk  r~>Rk 

Note  that  if  {Rk}  is  bounded,  then  there  exists  a  convergent  subsequence f?*  — >  R  .  We  will 
show  R  =  go  ,  implying  is  unbounded.  In  doing  so,  we  will  prove  the  existence  of  entire 
large  solutions.  Notice  that  since  vk  ,  uk ,  and  H  from  (4.3)  are  non-decreasing,  we  have 
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t  L 

vk  (r)  =bk+l  ^ "  J  s"1q(s)g(uk  ( s))dsdt 
0  0 
r  t 

<bk+^  tx~n |  sn~xq(s)H{uk  ( s))dsdt 
0  0 

oo  t 

\  bk  +  H(uk  (r)) J  tl~”  | sn  ]q(s)dsdt 


<i 


< Li. 

a,. 


0  0 


r  i  w  i 

ak  +  J  t'  iV | sn~lp(s) f  (vk ( s))dsdt  +  H(uk  (/*)) J  tl~n  J  s"  lq(s)dsdt 


V  o  o 


o  o 


-  Cxuk  (r)  +  C2H (uk(r))). 


The  value  C,  is  any  bound  on  bk  /  ak  ,  which  exists  since  a  ^  0 ,  and 


(4.33) 


CO  t  -.00 

C2  =  J  tx  "  ^ sn~x q{s)dsdt  <  - — —  ^sq{s)ds 


0  0 


<  oo . 


This  inequality  for  C2  is  obtained  by  integrating  as  in  (3.8)  and  (3.10)  and  is  finite  from  (2.10). 
We  then  define  h(t)  =  H(Cxt +  C2H(t)) .  Using  the  results  of  Lemma  4-1,  we  have 


h(  0)  =  H(  0  +  C2H{  0))  =  H  (0)  =  0, 
h(s )  =  H{Cxs  +  C2H(s ))  >0,  s  >  0, 


(4.34) 


and,  choosing  Cx>\, 


oo  |  t 


i  Lo 


-1/2 


w  i 

dt  =  J  J//(c15+c2//(5)) 

l  Lo 

00  t 

A  jffw 


-1/2 


dt 


i  lo 


dt 


<  oo. 


(4.35) 
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That  is,  h  satisfies  (2.12)  and  the  Keller-Osserman  condition  (2.2).  Lemma  1  in  [12]  then 
guarantees  that 


(4.36) 


is  well  defined  for  all  s  >  0 .  Note  that 


F'(s)  =  —  <  0  and  F"(s)  =  h^S\  >  0  . 
Ks)  [/?(*)]- 


(4.37) 


We  have  from  (4.33)  that 


A uk  =p(r)f(vk ) 

<  p(r)H(vk) 

<  p{r)H  [C{uk  (r)  +  C2H(uk  (r))] 
=  p(r)h(uk). 


Then,  using  (4.37),  we  calculate 


AF(uk)  =  F'(uk)Auk+F"(uk)\Vuk\ 
-1  h'(uk)  |  ,2 

\h(uk)\ 

p(r)h(uk) 


> 


h(uk ) 

-1 


h(uk) 
=  ~p(r). 


Rewriting  the  Laplacian  in  radial  fonn  and  multiplying  each  side  by  rn  1 ,  we  obtain 


( 

V 


r 


n  -1 


v 


>  -rn  1  p(r)  . 
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Integrating  over  [0,  r ]  where  0  <  r  <  Rk  gives  us 


d_ 

dr 


F{uk)>-r 


i 


s”  1  p{s)ds . 


Next,  we  integrate  over  [ r ,  Rk  ] .  Notice  that  since  uk  ( r )  — >  oo  as  r  — »  Rk ,  we  see  from  (4.36) 
F(uk  (r))  — »  0  as  r  — >  Rk .  Thus  integration  yields 


Kk  r 

- F(uk(r ))  >  -  J  tn~l  j sn~l p(s)dsdt . 

r  o 

That  is, 

Rk  r 

F(uk(r ))  <  J  tn~x  lp(s)dsdt . 

r  0 


Since  F'(s)  <0  for  s  >  0 ,  we  have 


uk(r)>F 


Kk 

|  tn-[\sn~'p(s)dsdt 

V  '■  0  j 


Now,  we  let  ^  >  co  so  Rk  —>  R  and  uk  —>  u  .  We  have 


I 

F-1  | rl\sn-lp(s)dsdt 


.  r  0 


<  w(r) . 


Letting  r  — >  R ,  and  since  limF(s)  =  0  implies  lim  F  '(s)  =  co  ,  we  have 

s-*co  s->0+ 


limF-1 

r->R 


R  r  X 

J  t"~l  jV' lp(s)dsdt 


J 


=  lim  F  1  (5)  =  00  <  lim  u(r) . 

i->  0+  r-*R 
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However,  recall  u{\  x  |)  and  v(|  jc  |)  have  central  values  (a,b)  eE  cz  S  and  are  entire.  This 
implies  u(r)  exists  for  all  r  e  [0,go)  .  Hence  R  =  qo  ,  and  our  proof  is  complete.  □ 

4.3  Examples 

In  this  section  we  present  several  examples  of  systems  that  satisfy  the  restrictions  given 
in  our  analysis.  Systems  of  the  form  (1.2)  which  satisfy  the  hypotheses  of  our  theorems  and 
lemmas  include  equations  such  as 


| A u  =  e  ^(vev),  |a»  =  2  3^v2, 

[Av  =  e^(ue“),  [Av  =  4~5^«3, 

In  each  of  these  examples,  the  functions  /  and  g  are  monotonic.  In  our  results,  we  also  have 
shown  existence  of  entire  large  solutions  for  systems  containing  non-monotonic  functions.  For 
example, 


|Au  =  2“H(4ev-4), 
{Av  =  3_W(5e“-5). 


A u=e  ^  {lv  +  |4vsin(v)|-l), 

< 

Av  =  e  ^  (2“  +|4wsin(«)|-l), 

also  satisfy  the  hypotheses  of  our  theorems  and  lemmas.  The  most  difficult  condition  to  check  is 
that  /  and  g  satisfy  the  Keller-Osserman  condition  (2.2).  We  guarantee  the  examples  in  (4.39) 
satisfy  this  requirement  since 

f(s)  =  g(s)  =  2  s  +|4jsin(j)|-l  >  2s  -1, 

/ (s)  =  g(s)  =  2s  Isin5|  +  (2S  -1)|cosj|  >  2s  -1, 


A u  =  3  'fy2v|sinv|  +  (2v-l)|cosv|j, 
Av  =  4^  {lu  |sinw|  +  (2“  -l)|cosw|), 


(4.39) 
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and  2'  - 1  satisfies  (2.2).  These  are  only  a  few  examples.  If  we  define  our  systems  piecewise, 
we  can  create  a  wide  variety  of  complicated  systems. 

Next,  we  shall  consider  several  of  these  examples  numerically.  Consider  the  simple 

system 


| A u  =  e  ^(vev), 
|av  =  e^x\ueu). 


(4.40) 


in  M5 .  From  our  analysis,  we  know  there  exists  a  set  of  central  values  S  for  which  entire 
solutions  exist.  Also,  we  know  these  solutions  are  large  for  central  values  in  E  given  in  (4.32). 
We  will  numerically  solve  this  radial  system  by  applying  a  Runge-Kutta  algorithm.  We  first 
solve  (4.40)  for  central  values  (0, 0)  e  S  .  Then,  holding  v(0)  =  0  constant,  we  steadily  increase 
the  central  value  u{ 0) ,  and  numerically  solve  for  each  iteration.  Figure  4-2  shows  these 
solutions  for  u{ 0)  ranging  from  0  to  2.  The  dashed  lines  represents  u  while  the  solid  lines  are 
solutions  v .  Figure  4-3  is  a  closer  examination  of  the  transition  from  entire  bounded  solutions  to 
large  solutions  on  bounded  domains.  In  other  words,  Figure  4-3  shows  the  behavior  of  solutions 
as  we  cross  over  dS . 
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Figure  4-2.  Numerical  Solutions  of  (4.40)  for  Varying  Central  Values 


Figure  4-3.  Numerical  Solutions  of  (4.40)  for  Central  Values  Near  Boundary 
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The  solutions  for  small  central  values  appear  bounded  and  entire.  As  we  increase  the 
central  value  for  u  ,  our  solutions  increase.  Once  the  central  value  is  large  enough,  our  solutions 
blow  up  and  exist  only  on  finite  domains.  This,  of  course,  agrees  with  the  results  from  our 
analysis.  We  can  perfonn  similar  numerical  trials  and  vary  both  u( 0)  and  v(0) .  From  this  data, 
we  construct  the  set  of  central  values  S  for  (4.40)  as  shown  in  Figure  4-4. 


Figure  4-4.  Entire  Solution  Existence  Region  for  Central  Values  of  (4.40) 

We  have  also  numerically  calculated  u{ 0)  and  77(0)  as  given  in  (4.10)  and  (4.16)  allowing  us  to 
construct  Tx  and  T2  from  Corollary  4-2.  We  do  this  to  allow  comparison  to  Figure  4-1.  Clearly 
our  numerical  results  shown  in  Figure  4-4  agree  with  our  analysis  summarized  in  Figure  4-1. 

In  the  previous  example,  S  had  a  simple  shape  due  to  the  monotonicity  of  our  functions 
and  the  symmetry  of  our  equations.  We  will  also  numerically  examine  a  simple  system  in  which 
/  and  g  are  non-mono  tonic.  Consider  the  system 
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Figure  4-6.  Numerical  Solutions  of  (4.41)  for  Central  Values  Near  Boundary 


Our  results  appear  to  be  similar  to  those  obtained  in  our  first  numerical  example. 
However,  if  we  examine  these  solutions  closer,  we  discover  several  significant  differences.  For 
example,  our  sequences  of  functions  are  not  increasing  as  we  increase  our  central  values. 
Indeed,  see  Figure  4-7. 
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Figure  4-7.  Numerical  Solutions  of  (4.41)  for  Varying  Central  Values 

Also,  our  set  of  central  values  S ,  shown  in  Figure  4-8,  appears  quite  different. 
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Figure  4-8.  Existence  Region  for  Central  Values  of  (4.41) 
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The  dashed  line  connects  the  largest  central  values  of  the  form  (0,v(0))  e  S  and  (w(0),0)  e  S  . 


For  this  non-monotone  example,  S  does  not  appear  to  be  convex.  If  it  were,  the  dashed  line 
should  lie  entirely  in  S .  Indeed,  we  have  (0,2. 29670),  (1.19673,0)  e  S  ,  and  if  S  were  convex, 
we  can  choose  n(0)  =  0.2  and  calculate 

v  = -2229670  2 29670  =  1.91286 

1.19673 

so  that  (0.2,1.91286)  e  S  .  However,  if  we  solve  numerically  for  these  three  pairs  of  central 
values,  we  obtain  Figure  4-9. 


0  10  20  30  40  50  60  70  80 

|x| 

Figure  4-9.  Solutions  Showing  Set  of  Central  Values  for  (4.41)  is  Not  Convex 

Our  solution  for  (w(0),v(0))  =  (0.2,1.91286)  is  not  entire  while  the  remaining  central  values 
yield  entire  bounded  solutions.  Thus,  S  appears  nonconvex  according  to  our  numerical 


4-35 


calculations.  This  is  very  interesting  since  S  must  be  convex  for  systems  with  monotone, 
convex  functions  /  and  g  (see  [13]  and  [22]). 
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V.  Conclusion 


5.1  Summary 

Our  research  began  by  examining  a  few  open  cases  left  from  the  work  of  Smith  [26],  We 
examined  semilinear  elliptic  equations  of  the  form 

A  u  -  p{x)ua  +  q(x)up ,  reQcii"  (5-1) 

for  the  sublinear  case  0  <  a  <  ft  <  1 .  We  proved  in  Theorem  3-1  that  (5.1)  has  entire  bounded 
solutions  when  nonnegative  functions  p  and  q  are  locally  Holder  continuous  and  (3.1)  holds. 

We  also  showed  in  Theorem  3-2  that  (3.14)  is  necessary  and  sufficient  for  (5.1)  to  have  an  entire 
large  radial  solution.  These  proofs  are  extensions  of  Theorem  22  and  Theorem  23  in  [26], 

Next,  we  presented  several  examples  of  equations  which  satisfied  the  hypotheses  of  our 
proofs.  We  found  numerical  solutions  to  a  few  of  these  problems,  and  visually  demonstrated 
results  of  our  analysis. 

After  examining  single  equations,  we  considered  semilinear  elliptic  systems.  We  studied 
the  radial  problem 


\Au  =  p(jx  |)/(v), 

\  ,  .  .  (5.2) 

[Av  =  q(\x  |)g(w), 

Lair  and  Wood  [13]  considered  (5.2)  for  f(v)  =  va  and  g(u)-  u,! .  Cirstea  and  Radulescu  [6] 
and  Peng  and  Song  [22]  generalized  their  results  by  considering  /  and  g  non-decreasing.  We 
built  on  the  foundations  of  these  works  to  establish  the  existence  of  entire  large  positive  solutions 
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of  (5.2)  for  non-monotonic  /  and  g .  We  only  required  that  G  ,  given  in  (4.3),  satisfy  the 
Keller-Osserman  condition  (2.2).  In  our  analysis,  we  established  several  other  notable  facts. 

First,  we  demonstrated  in  Theorem  4-1  that  system  (5.2)  has  an  infinite  number  of  entire 
nonnegative  solutions.  In  a  corollary,  we  showed  that  entire  solutions  for  this  system  exist  even 
when  G  fails  to  satisfy  (2.2).  We  only  need  that  f  ox  g  satisfy  the  Keller-Osserman  condition. 

Once  we  established  existence  of  entire  solutions,  we  characterized  the  set  of  central 
values,  S  .  We  proved  in  Theorem  4-2  that  this  set  is  closed  and  bounded,  and  in  a  corollary,  we 
provided  a  rough  geometric  description  of  the  set.  We  then  considered  solutions  with  central 
values  outside  S 

Lemma  4-3  showed  for  any  central  values  (a,b)<£S ,  a  >  0,b  >  0  ,  a  solution  to  system 
(5.2)  exists  on  a  bounded  domain.  This  domain  is  dependent  on  the  central  value  chosen.  Then 
in  Lemma  4-4,  we  proved  that  on  the  largest  domain  for  which  these  solutions  exist,  the  solution 
is  large.  That  is  lim  u(r )  =  go  =  lim  v(r )  where  Rc  d  is  defined  in  (4.25).  We  then  considered 

r  ~*Rc,d  r  ~^Rc,d 

solutions  to  (1.2)  where  the  central  values  a,b>  0  lie  on  the  boundary  dS . 

In  Theorem  4-3  we  used  many  of  our  earlier  results  to  finally  establish  the  existence  of 
entire  large  solutions  to  system  (1.2).  We  showed  these  solutions  exist  for  central  values  in  E , 
defined  in  (4.32).  Our  more  significant  results  for  systems  are  summarized  in  the  following 
table. 
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Table  5-1.  Existence  of  Solutions  for  System  (1.2) 


A  u  =  p(\x\)f(v) 
Av  =  q{\x\)g{u) 


Entire  Domain 

0  =  M" 

Solution  Exists/ 

Yes 

Requirements 

b,c,d,e 

Large  Solution  Exists/ 

Yes 

Requirements 

a,c,d,f 

a)  G(s )  =  min{  f  ( t ),  g(t)}  satisfies  Keller-0 ssennan  condition  (see  (2.2)). 

s<t 

b)  f  ox  g  satisfies  Keller-0 ssennan  condition  (see  (2.2)). 

c)  / (0)  =  g(0)  -  0 ,  f(s)  >  0  and  g(s)  >0  for  s  >  0  . 

OO  00 

d)  Functions  p,  q  not  identically  zero  at  infinity  and  J  tp(t)  <  oo,  J  tq(t)  <  oo 

0  0 

e)  Infinitely  many  solutions.  Central  values  form  closed  and  bounded  set  S 
(see  (4.1)). 

f)  Infinitely  many  solutions.  Central  values  on  E  (see  (4.32)) 


After  our  analysis,  we  identified  several  example  systems  that  satisfied  the  requirements 
of  our  theorems  and  lemmas.  We  then  examined  numerical  solutions  for  several  of  these 
examples.  In  addition  to  demonstrating  the  results  of  our  analysis,  our  numerical  results  suggest 
the  set  of  central  values  is  not  convex. 


5.2  Further  Work 

While  single  equations  similar  to  (1.1)  have  been  studied  extensively,  there  are  still  many 
open  problems.  Some  of  these  potential  areas  of  study  can  be  seen  from  Table  2-1.  Condition 
(c)  from  this  table  has  been  shown  to  be  necessary  for  the  existence  of  solutions.  It  is  unknown 
if  other  conditions  may  be  relaxed.  For  example,  it  is  unknown  if  entire  large  solutions  exist  for 
the  superlinear/mixed  case  when  condition  (a)  does  not  hold.  That  is  when 


r  max  p(x)dr  =  oo  or  r  max  q(x)dr  =  oo. 

J  |xl=r  »  |xl=r 
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In  fact,  a  similar  result  remains  unknown  for  the  single  term  equation  A u  =  p(x)ua  .  Also,  the 
existence  of  entire  large  solutions  for  the  sublinear  case  remains  unknown  if  condition  (b)  of 
Table  2-1  is  not  true,  that  is  when  (1 . 1)  is  not  radial.  We  believe  that  each  of  these  problems  will 
be  extremely  challenging  to  solve. 

A  much  simpler  problem  may  be  to  examine  properties  of  the  2-term  equation  and  extend 
them  to  an  arbitrary  n  -tenn  equation.  Results  will  likely  be  similar  to  the  2-term  problem,  but 
this  would  be  an  interesting  exercise. 

There  are  fewer  results  for  system  (1.2)  as  compared  to  the  single  equation,  and  therefore 
more  opportunities  exist  for  study.  Allowing  /  and  g  to  be  non-monotonic  has  created  unique 
problems  that  would  be  worth  examining  further.  For  example,  we  did  not  characterize  our  set 
of  central  values  beyond  closed  and  bounded.  For  monotonic  functions  /  and  g ,  [6]  and  [22] 

both  show  (a0  ,b0)eS  implies  ( a,b )  e  S  where  0  <  a  <  a0  and  0  <b<b0.  We  do  not  know  if  a 
similar  result  holds  for  non-monotonic  functions.  We  have  not  encountered  numerical  examples 
that  would  suggest  otherwise.  However,  we  considered  only  simple  systems  in  our  numerical 
trials.  Also,  we  are  unsure  if  large  solutions  exist  if  G  ,  as  in  (4.3),  does  not  satisfy  the  Keller- 
Osserman  condition  (2.2)  Indeed,  we  have  shown  entire  solutions  may  still  exist,  but  can  any  of 
these  solutions  be  large?  Finally,  we  would  like  to  further  relax  our  hypotheses.  Instead  of 
requiring  that  G  satisfies  the  Keller-Ossennan  condition,  we  would  like  to  require  that  only  / 
and  g  need  to  satisfy  this  condition.  When  we  began  studying  systems,  this  was  our  original 
hypothesis.  However,  to  use  our  arguments,  we  must  have  a  mono  tonic  function  below  /  and 
g  that  has  all  the  same  properties  as  G  .  Most  importantly,  we  need  the  function  to  satisfy  the 
Keller-Ossennan  condition.  At  this  time,  we  believe  such  a  construction  should  be  possible,  but 
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we  have  encountered  difficulties  when  trying  to  build  such  a  function.  It  would  be  worthwhile  to 
examine  these  areas  more  closely. 

These  are  just  a  few  examples  for  future  study.  In  general,  this  field  provides  many  open 
problems,  both  challenging  and  simple,  making  it  well  suited  for  students  and  researchers  of  any 
level.  Further,  the  wide  range  of  application  for  elliptic  theory  allows  this  area  of  study  to  be 
valuable  in  a  variety  of  fields. 
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